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Lubrication is a part of tribology/ the science and 
technology of interacting surfaces in relative motion and 
of related subjects and practices ♦ In the process of 
interaction/ part of the energy supplied to it is lost due 
to undesirable frictional forces opposing the motion. 
Lubrication_practice is an art of controlling these frictional 
forces/ reducing wear / distortion etc. by introducing some 
kind of lubricant in the separation gap formed by the 
interacting surfaces (bearing) • This lubricant plays a vital 
role during the process of lubrication of bearings* Bearings 
are very important for the successful development of the 
modern science and engineering. It seems rather impossible 
to imagine the progress of the advanced technology without 
the application of bearings* 



Separation between the bearing surfaces under the 
principle of hydrodynamic lubrication is achieved due to 
viscous shearing action of the fluid- This mechanism 
works very well in case of lubricated rolling bearings# 
which form an integral part of majority of machines* 

The lubricant# under the hydrodynamic lubrication of 
heavily loaded rolling bearing# is subjected to very high 
shear stresses and velocity gradients* The resulting 
extremes of pressure and temperature act on the lubricant 
for a very short duration which may cause the lubricant to 
behave as a non-Newtonian fluid [if 2 ]* In addition increasing 
use of additives (high molecular weight polymers) in the 
lubricant is also responsible for the deviation of its 
behaviour from Newtonian [ 3 ]* A number of fluid rheological 
models have been postulated to account for the non-Newtonian 
behaviour. One such model is that of the power law which 
has found extensive applications in recent years [ 4 ]« 

The lubricant temperature field is the consequence of 
hydrodynamic lubrication process which consists essentially 
of two surfaces in relative motion shearing a thin layer of 
viscous fluid film* This temperature field becomes stronger 
in case of heavily loaded contacts* Since the lubricant 
viscosity is strongly dependent on temperature# it becomes 
necessary to acquire knowledge of temperature field in the 
lubricant [ 5 ]* 



There are several other factors which can influence 
the temperature and pressure field in the lubricant* The 
fluid inertia is one of than which is generally small 
and may be neglected in comparison to viscous force- 
However* with the current trend towards design of bearings 
operating with high speeds* and use of lubricants having 
low viscosities and high densities* the lubricant inertia 
may attain significance for some range of moderately large 
Reynolds numbers [ 6 ]. This effect may become stronger 
in case of thick film bearing under squeezing motions [7]» 

Another important factor affecting the bearing lubrication 
characteristics is the compressibility effect* The compress- 
ibility effect may not be neglected for high speed flows 
or flows with large variation of pressure and temperature [ 8 ]. 
In fact the density is also sensitive to exceedingly high 
temperatures and pressures* 

The classical lubrication theory constitutes the 
solution to the standard form of Reynolds equation. However* 
much of the modem lubrication theory deals with the departure 
from it* This departure may be due to the thermal effects* 
inertia forces* compressibility* undetermined film thickness 
and non-Newtonian behaviour of the lubricant (as mentioned 
earlier) * The present thesis work is devoted to the study of 
some of these important factors* The thesis consists of a 
general introduction and five chapters- 



General introduction is introductory in nature* It starts 
from the field of terminology followed by history of lubricated 
bearings and a literature survey (theoretical and experimental) 
of lubricant rheology/ fluid rheological models/ thermal 
and inertia effects and compressibility. It may be noted 
that attention has been focussed mainly on the problem of 
lubricated rolling/sliding contacts* 

In chapter 1/ a theoretical problem of pure rolling of two 
identical rigid cylindrical rollers moving with equal speed 
and lubricated with an incompressible power law fluid is 
considered under the usual assumptions of hydrodynamic 
lubrication. The lubricant consistency m is assumed to vary 
exponentially with pressure and temperature* The moving 
boundaries are assumed to be adiabatic* The modified 
Reynolds and energy equations are derived and are solved 
numerically for various values of the flow index n and the 
squeezing parameter q* Pressure/ temperature and the film 
cavitation point are computed* It is observed that the 
effect of temperature on consistency is to shift the position 
of pressure peak slightly away from the centre line of contact- 
Pressure and temperature increase significantly with n/ 
however both decrease as q increases* The consistency change 
in the pressure peak region is quite significant. 
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Chapter II deals with a similar geometry followed by a 
similar theoretical approach* The difference appears here 
due to the fact that the lubricant is assumed to be 
compressible* The lubricant density is assumed to vary 
with pressure as well as temperature- In view of the 
extra complication which arises due to the consideration 
of compressibility effect/ squeezing effect has been 
dropped* Pressure / temperature/ load and traction etc* 
are calculated for compressible as well as incompressible 
fluids- Compressible load and traction are found to be 
comparatively higher than those for incompressible case- 
The compressibility effect is found to be significant for 
Newtonian as well as dilatant fluids* 

In chapter III also/ the geometry of the problem and 
the mathematical formulation remain the same/ as in chapter 1/ 
however/ the adiabatic assumption is relaxed here* The 
heat of convection which is comparatively smaller to that 
of conduction is neglected* The lubricant consistency is 
assumed to be dependent on pressure and the mean film , 
temperature* Various bearing characteristics such as pressure 
the mean film temperature/ load and traction are calculated 
as functions of n and q. The bearing characteristics depict 
a similar trend with respect to n and q as in chapter I* 



The problem dona in chapter IV is similar to that of 
chapter III. The difference here is that the fluid inertia 
and convection terms are retained in the momentum and the 
energy equations respectively. The presence of inertia 
and convection terms makes the solution very complicated. 
Thus averaged values of both the terms have been considered. 
The results for various bearing characteristics with and 
without inertia have been calculated. It is observed 
that the inclusion of inertia force leads to increase in 
pressure and thus also the load. Squeezing motion is found 
to have significant effect on inertia* 

The last problem done in chapter V is also a generalized 
version of the problem formulated in chapter III. The major 
difference here is that both the rollers are assumed to have 
different dimensions and velocities/ squeezing effect is 
here dropped . Various important results are calculated 
for various values of n and sliding parameter U- 
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NOMENCLATURE 


(n/(2n+l)) n 

1/«P S 

n 2 /( 2n+l) ( 3 n+1 ) 

2(2n+l)/(3n+2) 

specific heat of lubricant at constant volume 

((4n+2) /n) n (U/h ) n (2R/h J 1 ' 72 
' ' o ' o 

specific heat of the lubricant at constant pressur 
n./((?n+l)2 n+1 ) 

a number as defined in (4«15) 

2 

dimensionless number ( = U p/c) 

X 2 -xJ+2g(X-!-X ;L ) (Chapter IV) 

H-P*h' 3J '/P 1 (Chapter II) 

X 2 -X 2 +2q(X+X 1 ) (Chapter I) 

H-P H* e -APi+e(Ti~T 0 ) (chapter II) 

X 2 -X 2 +2q(X+X 1 ) 

-F 


a number as defined in (4*15) 

-h+p*h*/p 2 

-H + P^H%"^ p 2 +e ^2“^o^ 
film thickness 

minimum film thickness 

film thickness at the point of maximum pressure 

film thickness at x = x 2 

l+X^ 

1-tX* 2 
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1+X 1 

inertia term defined by relation (4*23) 
viscous term as defined in (4-2?) 
viscous term as defined in (4*33) 
lubricant thermal conductivity 

*| /O 

(Uh o /(4ka) ) (h o /2R) ' 
lubricant consistency 

consistency at ambient pressure and temperature 
consi stency ( =2mc n « ) 

consistency at ambient pressure and temperature 

( =2m c a ) 
o n 

flow behaviour index 
hydrodynamic pressure 
pressure! =«p) 

bulk modulus of compressibility 

(ch 0 /4kau) (h Q /2R) x/2 

squeezing parameter (= ^(R/2h Q ) ^ 2 ) 

radius of the cylinder 

equivalent radius 

l+X 2 -2q(X+X 1 ) 

lubricant temperature 

temperature ( ={? (T-T o ) ) (Chapter i) 

temperature (=j3T) (Chapter II* IV and V) 

traction force 
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dimensionless traction force (=-T f /(h /2^) ) 
surface temperature (constant) 
temperature (=T^/k^) (Chapter III) 

(*0T h ) (Chapters IV and V) 

ambient temperature 

temperature ( -T 0 /k^ ) (Chapter III) 

temperature (~{3T o ) (Chapter IV and V) 

surface temperature at y = o 

* T ol 

the mean film temperature as defined in relation (3*6) 

mean temperature (=T m ,/k^) (Chapter III) 

mean temperature ( =0' T m ) (Chapter IV and V) 

lubricant viscosity in x-direction 

u/U 

the mean velocity as defined in (4»14) 
velocity of the cylinder in pure rolling 

U l- /tJ 2 

velocity of the plane surface 

velocity of the cylinder under rolling/sliding 
contacts 

velocity of the lubricant in y-direction 
velocity <, =( 2v/U ) ( 2R/h Q ) ) 

V/2 is the normal velocity of the cylinder 
load in x-direction 
dimensionless load (=W /(2h AO) 

-A. O 



1 /? 

dimensionless load (=W /((2Rh ) - V")) 

Y' O 

co-ordinate axes 

-x* or -x^ is the point of maximum pressure 
cavitation point 

dimensionless distance (=x/(2Rh ^Z 2 ) 

• o 

dimensionless distance (=2y/h ) 

1+X* 2 

1+X 2 

cavitation point (=1+X~) 
pressure coefficient 
temperature coefficient 

P r E t (n - /(3n+l) ^ 
p/(Pccc) (Chapter i) 

P /(P o c p a) (Chapter II) 

PJE. (n/(4n+l)) 
r t: 

U 2 h o^ <V 2R)1/2 / CMx) 

locus of points at which lubricant velocity 
gradients vanish and is measured from y = o 
6.^0 

a value of & at X = 
values of 6 at x = +x^ 

**/h 0 

pu 2 a 

lubricant density 
lubricant density at x = -x* 
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P Q the lubricant density at ambient pressure and 

temperature 

p p /p c 

Pj density of the lubricant X = 

e temperature coefficient 

e e/p 

Subscripts • 

1.2 Usually refer to the respective quantities in' 

the inlet and the outlet regions separated by the 
point at which velocity gradient vanishes otherwise 
explicitely mentioned* 

Bar denotes dimensionless number* 
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GENERAL INTRODUCTION 


1* LUBRICATION OF BEARING, /FIELD TERMINOLOGY I 

Lubrication is a part of ' tribology* / the science and 
technology of interacting surfaces in relative motion and 
of related subjects and practices [l]» As the interaction 
between surfaces occurs during the process of sliding/ rolling 
etc* / part of the energy supplied to it is lost due to 
undesirable frictional forces opposing the motion* Lubrication 
practice is an art of controlling these forces/ reducing 
wear/ distortion etc*/ by introducing some kind of substance 
between these interacting surfaces having a relative motion* 
Insertion of such substance between the mating surfaces makes 
them operate more efficiently/ i. it serves to separate 
rubbing surfaces/ aids in the removal of heat/ prevents 
rust formation/ protects bearing from dirt and minimizes 
maintenance etc*. This substance is called a lubricant* 

The most commonly used lubricants are oils and greases* 

However/ it is not unusual to use other liquid substances/ 
as well as solid and air as lubricants* 

Bearing is a machine component or a part whose function 
is to support a second member / preventing its motion in the 
direction of applied load but at the same time allowing 
motion in another predetermined direction [2]» It is not a 
matter of surprise that the longstanding appreciation of the 
importance and applications of bearings are quite unlimited* 



2 


It seems impossible to imagine the progress of. modem 
science and technology without the application of bearing 
which forms an integral part of the majority of machines* 

One look at the devices and machines which typify present 
day civilization is enough for realizing how important 
successful bearing operation is to our daily living* 

The general theme which runs through the consideration 
of bearings of whatever type is that a bearing must be 
considered from three aspects t no one of which is completely 
separable from others* These aspects are * (l) the design 
of bearing \ ( 2 ) the material of bearing / and ( 2 ) the 
lubricant (used)* No bearing is completely designed without 
proper consideration and integration of these three aspects [ 3 ] • 
Bearing analysis (especially rolling bearing) together with 
the third aspect will be the point of discussion within 
the frame work of this thesis* 

The rolling contact bea£p.ngs utilize the rolling action 
of balls and/or rollers having point or line contacts / and 
hence exercise minimum rolling frictional resistance. 

Bolling contact bearings are usually employed to allow 
either rotary motion or relative linear motion and sometimes 
a combination of both the types of motions [4]- A historical 
development has been given in the subsequent section. 



2* HISTORY OF LUBRICATED BEARING AND RECENT DEVELOPMENTS *. 

The art of lubrication was known many centuries ago 
and lubricants such as animal fats / vegetable oils and even 
mineral oils were available* The chariot wheel* axle of 
which had been lubricated with a sticky greasy type substance 
( during the era of Egyptian Pharaohs) / perhaps reminds of 
the early days of lubrication [5j* It was not until the 
end of 19th century that the true nature of lubrication 
mechanism was fully understood and a new branch of fluid 
mechanics / known as lubrication mechanics * was formed [6j. 

Him [7] who may be called the father of lubrication* discovered 
film lubrication in 1854* But the first research credit goes 
to Tower in England and Petroff in Russia [8]» An experiment 
on oil lubricated bearing conducted by Tower [ 9] in reply to 
a mistaken idea provided a break through in this field* He 
explained that the interacting surfaces might be separated 
by an oil film without metallic contact. He further concluded 
that a loaded oil bearing might be subjected to local pressure 
substantially higher than its mean pressure- This experiment 
attracted the attention of Prof* Reynolds [lo] who was the 
first to give a complete mathematical formulation related 
to Tower* s experiment* This was followed by Sommerfeld [ll]« 
Ever since * this field has never lost its momentum because 
of its wide application in science and technology. Since 
the onset of the atomic age / bearing lubrication and other 
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scientific progress has occured virtually at an exponential 
pace* Today in the space age / continually increasing 
demands are being made on various sophisticated engineering 
equipments containing bearings* This has necessitated a 
reappraisal of the bearing design procedure in order to 
minimize frictional resistance and mitigate other undesirable 
phenomena* 

Rolling element bearing is of relatively recent origin/ 
so far as its common engineering application is concerned/ 
dating back roughly to the year 1900 [3j» The principles 
involved were known many centuries ago but it is only with 
the advent of modern materials and techniques/ the rolling 
element bearing has found extensive use* Since 19 60/ 
substantial attention has been given to the mechanism of 
rolling bearing lubrication and rhe rheology of lubricants* 
Manufacturers for world wide markets increased substantially 
in the seventies which have served to provide standard 
bearings of outstanding performance an relatively low cost 
to the consumers [ 4 ]* The advan cages/ disadvantages over other 
types of bearings are clear because rolling friction is far 
less than that of others* This has been discussed by Wilcock 
and Booser [ 3 ] * 

Mathematical treatment, of lubricated bearing is quite 
complicated* In the early days researchers working in this 
field were acquainted with the fact that the assumptions/ 
such as no slip/ isoviscous/ isothermal/ Newtonian lubricant 
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etc*/ are not valid for bearing under severe operating 
conditions* Consequently their theoretical results were 
not in good agreement with the experimental observations* 

It was thus realized that those assumptions cannot be made 
in order to get closer agreement with the experiment * 

With the advent of fast digital computers and various 
numerical techniques (FDM# FEM/ BEM etc*)/ mathematical 
treatment became easier and analysts began to relax assumptions 
and more realistic solutions were obtained* However even 
today there exists no mathematical treatment which takes 
into account all the factors likely to influence lubrication 
characteristics* A theoretician must therefore resort to a 
gradual relaxation in the assumptions- 

3* LUBRICANT RHEOLOGY * 

Rheology is the science dealing with the deformation 
and flow of materials* Its significance in the field of 
lubrication has arisen in part from the realization that 
the lubricant/ especially between gears and rollers/ is 
subjected to quite high shear stresses and velocity gradients* 
The rheological properties of lubricant are therefore likely 
to have important effects on its actions [12]* Serious 
attention was paid to the rheology of lubricants in the 
sixties however the mechanism of bearing lubrication is not 
yet fully understood. It is also known that rheological 
properties of the lubricant in the proximity of the solid 
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surface change significantly. Particularly this phenomenon 
is evident when the additives in the lubricant are surface 
active in the EHD regime [13-14]. 

With the increasing demands being placed on lubricants / 
for automotive engines/ transmission and high speed loaded machinery 
much attention has been paid to the development of improved 
lubricants (mineral oils with additives) « The first successful 
engine lubricant additive was produced by California Research 
Corporation in 1935 when the addition of Aluminium Naphthenate 
was successful in overcoming the sticking of Caterpilar diesel 
engine [15J. A small percentage of these chemical additives 
IS incorporated to satisfy those demands that the natural 
properties of mineral oils cannot adequately cope with* 

Additives are used to fulfil specific requirements such as 
anti-oxidants (Sodium soap)/ anti-wear (Calcium soap)/ E P 
additives (Lithium soap)/ corrosion inhibitors (Aluminium soap)/ 
friction modifiers (Barium soap) / metal daactivators 
(Aluminium complex) / pour— point depressants (Bentonite clay) / 
tackiness additives (silica) and water repellants (Carbon/ 

G raphite) etc . . 

Kingsbury [16] may be considered to be the pioneer as 
far as the rheological abnormalities in lubrication is 
concerned. He observed enhancement of viscosity of the fluid 
in the region of attraction of the surface molecules of the 
metals* Needs' [17J experiment on rheological behaviour is 
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of great importance and interest- He squeezed fluids 
between two steel plates and found that a stable residual 
film was formed. Though some of his observations are question- 
able/ it appears reasonable that an alteration of liquid film 
occurs near the bounding solid surface which makes it 
substantially more viscous compared to the bulk of the liquid- 
Fuks [18- 19] also conducted a number of experiments in 
support of N eeds # observations/ using an apparatus similar 
to that used by Needs- The works of Fuks and Needs were a 
matter of controversy and criticism/ mainly by Hayward and 
Isdale [20]. They pointed out that the asperities on the 
solid surface might account for residual film. The dirt 
and the undissolved material in the fluid film have been 
also offered as an explanation. More recently Cameron et al 
[21-24] studied the influence of surface active compounds in 
viscous lubricating systems* The lubricants were pure Paratfins/ 
mainly Hexadecane and other long chain polymer compounds* 

It was concluded that the oil film was influenced by the 
additive and solvent- Paul and his co-workers[ 25-27] and 
others [ 28 ] examined lubricant rheology experimentally on the 
basis of structural rearrangement under high shear stresses* 

It was convincingly demonstrated that the change in rheological 
properties comes out to be significant. Hirst and his co-workers 
[29-30] also examined lubricant rheology during impact 
theoretically and experimentally/ using mineral oils and 
silicon fluids. He came up with a conclusion that these fluids 
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exhibited non— Newtonian behaviour (non-linear relation 
between shear stress and the rate of strain) » Further/ 
using a model for non-Newtonian fluid/ he showed a good 
agreement with the experiment* Walter [3l] presented a 
survey on the development on theoretical and experimental 
rheology which provides an adequate justification for the 
increasino importance of the rheology in lubrication studies* 

4* FLUID KHLGLOGICAL MODELS * 

It is well known that many mineral oils which are most 
commonly used as lubricants manifest non-Newtonian behaviour 
under severe operating conditions* The changes in the 
stress and the rate of strain are so rapid that the viscous 
linear model (Newtonian) is no longer adequate [32-35]* 

Johnson [36] pointed out that if the shear stress in a. typical 
lubricant exceeds S?10^ >Pa/the fluid exhibits non-Newtonian 
behaviour. This behaviour is usually encountered under 
heavily loaded conditions* It was pointed out earlier 
that the additives (such as high molecular weight polymers) 
are used in mineral oils to obtain or enhance many properties 
that the ordinary mineral oils do not possess* The presence 
of these polymer additives in mineral oils is responsible 
for deviation from the Newtonian behaviour* Non-Newtonian 
behaviour may also bo displayed due to the presence of 
impurity in the lubricant* Further the dirt particles collected 
during the process of lubrication and manufacturing may make 
the operating lubricant non-Newtonian* 



9 


A number of empirical fluid models have been suggested 
to account for non-Newtonian effects in lubrication- A 
Bingham plastic flow model/ a first approximation to describe 
the behaviour of grease lubrication/ has been used extensively- 
Milne [37] presented a theoretical study of grease lubricated 
slider bearing using this model- His analysis showed the 
existence of cored regions in the grease film. Sasaki et al 
[38] also used the same model and presented a theoretical 
study of grease lubricated roller bearing on the basis 
of cored regions- Wada et al [39] studied EHD lubrication 
between two cylinders using the same model where the plastic 
velocity and yield stress of the lubricant were assumed to 
be dependent on pressure. Experimental work on velocity 
determination in a grease lubricated cylindrical sliding/ rolling 
contact has also been reported. Mutuli et al [ 40 ] made 
comparison between the theoretical results of a parabolic 
velocity distribution similar to that of a Newtonian fluid# 
derived from a Bingham plastic flow model. 

Burton [4l] studied the lubrication of rollers using 
Newtonian and non-Newtonian viscoelastic Maxwell model* A 
comparison of results (pressure distribution) for both the 
models indicates that shear elasticity can hav-e effects of 
first Order significance on both the magnitude and form of 
pressure distribution functions reducing load carrying capacity 
and maximum shear stress in the film by sizeable amounts- 
Walter [3l] modeled the polymer thickened oil as elastoviscous 
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fluid and made various observations on lubrication of gears 
and rollers* Johnson and Tevaarwerk [42] formulated a new 
rheological equation on the lubricant by addition of a linear 
elastic element to non-linear viscous element* This model 
was further extended to include -cnemial effects by Houpert [43j? 
Bair and Winer [44] offered a modified Maxwell model that 
yields good agreement with the experiment* Heyes and 
M.ontros« [45] also made a theoretical study of SHD small 
strain behaviour of lubricant using viscoelastic model that 
involves reformulation of the Maxwell model in terms of 
Volterra convolution integral equation- Phan-Thien [46] 
studied flow of viscoelastic fluid squeezed between two parallel 
plates* It has been shown that the difference between two 
forces acting on the top and bottom plates is entirely an 
inertia phenomenon and arises due to the acceleration of the 
top plates (bottom plate is stationary). He and his co-workers 
also studied the squeezing problems with Maxwell fluid [47]/ 
Oldroyd fluid [48]/ and real elastic liquid [49]* 

Hsu and Saibal [50] used the. constitutive relation (cubic) 
which is applicable to pseudoplastic and dilatant fluids* 

They found that for the pseudoplastic case the load capacity 
is less than that of Newtonian* Wada and Hayasi [5l] analysed 
the hydrodynamic lubrication of a journal bearing by pseudo- 
plastic fluid both theoretically as well as experimentally* 

They also concluded that the load capacity/ drag force and 
pressure distribution are reduced when compared to the 
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performance of Newtonian lubricant* Further Wada et al [52] 
applied the same model for cylindrical rollers including EHD 
effect and arrived at a similar conclusion. Very recently 
Wada and his co-workers [53] applied the same constitutive 
relation for a parallel circular squeezing film bearing 
including the affect of inertia and obtained various bearing 
characteristics considering sinusoidal squeezing motion* 

Among the various models postulated to describe the 
flow behaviour of a non-Newtonian fluid# the power law 
model has perhaps found the most extensive application in 
various fields [54-58]. It is also probably the most criticized 
and most maligned model in all of rheology [59]. The appeal 
of this model is however evident when n / the flow behaviour 
index * is unity the model reduces to the Newtonian. For 
0 < n < 1 * it describes a rheogram with the shape of a 
pseudoplastic fluid. For n > 1 # the curve characteristic 
of dilatant behaviour is found. Furthermore# from the values 
of n one can readily obtain a qualitative picture of the 
degree of shear- rate dependence shown by a given material* 

The breadth of the application is shown in table 1 where a 
list of power law constants has been compiled [59-61] * 
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Davson and Wilson [62] were perhaps among the earliest 
researchers who used the power law model in lubrication. 

They examined the EHL capabilities of silicon fluids (modeled by 
power lav/) and concluded that they are poor lubricants for 
heavily loaded sliding surfaces due in part to the thinner 
EHD films formed as a result of their non-Newtonian behaviour- 
Safar and Shavki |_63] studied a general solution for pressure 
distribution and load capacity for a hydrodynamic thrust bearing 
operating with power law lubricants* It was established 
that the load capacity of the bearing may be increased several 
folds with increasing value of the viscometric exponent n. 

Safar [64-65] applied the same model for journal bearing 
under steady and dynamic loading conditions and arrived 
at the same conclusion. Sinha at el [ 66] studied the lubricat- 
ion of rollers under lightly loaded condition using this power 
law model and analysed load and traction ratios as functions 
of flow index n and squeezing velocity ratio q. A table 
for the points of cavitation was also presented with respect 
to n and q* Further Sinha and Raj [67] considered this 
problem for heavily loaded case where the consistency variation 
m was also accounted for. Salem et al [63-69] applied the 
same model for spherical bearing. They also emphasized'on the 
consistency variation and extended the same problem for non— 
isothermal case* Dein and Elrod [70] formulated a generalized 
Reynolds equation for this model/ applicable to journal bearing* 
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Blkouh at al [7l] studied the laminar squeezing between 
two annular surfaces using this model* Further they generalized 
this to include the effect of difference in pressures at inner 
and outer boundaries [72]. A theoretical study of short 
journal bearings performance characteristics for the same 
fluid model was presented by Buckholz [73]* Fluid film pressure 
profiles in short bearing of arbitrary azimuthal length has 
been obtained using matched asymptotic expansions in the 
slenderness ratio* Further he and Hwang [74] discussed the 
accuracy of the short bearing approximation applicable to 
Newtonian as well as non-Newtonian power law lubricants* 

Raj [75] used this model for various geometries in his entire 
dissertation work. Nailwal |_7 6j uisb used this model to study 
seal behaviour. 

Apart from this/ there are several other non-Newtonian 
models which have been used theoretically and experimentally 
and it is beyond the scope of this thesis to refer to all the 
non-Newtonian models* 

5* THERMAL EFFECTS ‘ 

With the present day tendency towards limit design / it 
is becoming increasingly important that the bearing operating 
characteristics be predicted accurately* Because of the strong 
dependence of lubricant viscosiiy on temperature/ it becomes 
necessary to acquire knowledge of the temperature field in 
the lubricant [77]. 
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The temperature field in the lubricant is the consequence 
of hydrodynamic lubrication process which consists essentially 
of two surfaces in relative motion shearing a thin layer 
of viscous fluid film* With the exception of bearings 
running at low speeds/ the heat balance is an important 
part of the analysis of bearing operating with fluid film • 
The operating temperature of the bearing is also of interest 
because the fatigue strength of the bearing material is a 
function of temperature* 

In order to account for thermal effects/ the classical 
isothermal theory of hydrodynamic lubrication has to be 
extended. Energy balance equation is required and has to 
be solved simultaneously along with the conservation equations 
of mass and momentum. Coupling of these governing equations 
makes the solution very complicated and further simplification 
(approximation) may be necessary. 

Another approach for the consideration of thermal effects 
in lubrication has also been reported- Viscosity-temperature 
relation has been replaced by viscosity film thickness 
relation on the assumption that the highest temperature occurs 
in the zone whssre the film thickness is least [l4/78j . This 
approach does not incorporate energy equation at all/ so the 
complexity of the problem made by energy equation is overcome. 
But on the other hand/ the temperature field in the lubricant 
is not fully known. It seems therefore that this way of 
approach has been used sparingly* 
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A number of papers [79-82] touched upon the effects of 
temperature on hydrodynamic bearing performance and requirements* 
Tipei and Nica [ 83 - 84 ] obtained the three dimensional 
temperature variation of journal bearing oil film* Separate 
relationship were established for both the convergent and 
divergent regions of the bearing taking into account viscosity 
variations and side leakage* The theory was found to be in 
good agreement with the experimental data obtained* McCallion 
et al [85] studied the thermal effects by obtaining solutions 
to momentum and energy aquations for adiabatic and THD 
conditions* They compared the results with isothermal case 
and found that THD solution lies somewhere between isothermal 
and adiabatic. Crook [86-87] considered the effect of 
temperature on lubrication of rigid cylindrical rollers by a 
Newtonian incompressible fluid theoretically and experimentally* 
Pressure/ temperature distributions/ film thickness and traction 
curves were plotted for various parameters* Cheng and 
Stemlicht [88-89] considered very general case of non-iso thermal 
lubricating layer/ presenting an account of viscosity and 
density changes over the length and thickness of the film 
between two non-rigid cylinders rotating with unequal speeds * u 
Numerical solutions to the coupled Reynolds and energy equations 
were presented for hydrodynamic conditions* Dayson [90] studied 
the affect of non-Newtonian linear and non-linear Maxwell 
fluid model on the rolling/ sliding contact* He signified 
the region of contact in terms of linear/ non-linear and thermal 
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regions and presented the variation of apparent viscosity 
with pressure and rolling speeds and made comparison between 
linear and non-linear models* A case of EHD traction for a 
very particular non-Newtonian model was analysed by Kannel 
and Walowit' [91 J • They compared their results for traction 
with Crook's results (considering viscosity variation with 
pressure and temperature)- Johnson and Greenwood [92j 
extended Crook 7 s result on the effect of shear heating on 
EHD traction/ using an Eyring fluid rather than a Newtonian 
fluid* Conry [93] solved thermal EHD problem for some fluid 
to compare the predicted tractions with those obtained 
experimentally* Ghose and Hamrock [94] presented a numerical 
solution to the problem of thermal EHL of line contact using 
finite difference formulations* Pressure and temperature 
distributions and film shape for fully flooded conjunctions 
were obtained for a Paraffinic lubricant and various dimension- 
less speed parameters while the dimensionless load and material 
parameters were held constant. Houpert et al [95 ] analysed 
the rheological and thermal effects in lubricated elliptical 
Hertzian contacts with collinear speeds* Stress/ temperature 
and shear strain rate distributions aLong and across the film . 
were calculated for non-Newtonian visco-el as topi astic fluid 
including lubricant convection and conduction effects* Further 
Houpert [96] calculated traction force in EHD contact using 
Ree-Eyring fluid model and Roelands relationship [97] for 



estimating lubricant viscosity. The Roe-Eyring shear 
stress is found to be greater when the Roelands viscosity 
relationship is used- Szeri [98] in his most recent paper 
emphasized on THD effect in lubrication of gear teeth. He 
suggested that in bearing under severe operating conditions/ 
thermal effects may be significant to the extent that 
prediction of bearing performance is no longer possible when 
based on the assumption of a uniform film temperature- There 
are several other notable papers on thermal affects in lubrication 
for different geometries [99-101] and it is difficult to 
refer to all of them. 

6- COMPRESSIBILITY EFFECT * 

Lubricant compressibility is the measure of change of 
volume and density when the lubricant film is subjected to 
pressure and tension. This effect is caused by an abrupt 
change in pressure and temperature. Pai [102] pointed out 
that it cannot be neglected for high speed flows or flows 
with largo variation of pressure and temperature- It was 
supported by Dowson et al [l03j who also examined the 
influence of speed and lubricant compressibility on a theoretical 
solution to the EHD problem. Speed was shown to be the most 
important variable in the problem and the inclusion of the 
lubricant compressibility in the theory was found to have an 
appreciable influence on pressure distribution but a little 
effect on minimum film thickness* Since extremes of pressure 
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and temperature (together with large velocity) are often 
encountered in the lubrication of gear teeth and rolling 
contact/ this effect or phenomenon attains increasing importance 
in heavily loaded contacts* The effect of compressibility 
is overcome by a reduction in the film thickness which 
modifies the pressure and temperature quite considerably* 

Change in density might be expected to reduce the magnitude 
of the pressure recovery near the periphery of the bearing/ 
but accentuate phenomena near the inlet [ 104 ]* 

Houpert and Hamrock [105] who examined the compressibility 
effect of EHD lubricated line contact by Newtonian fluid 
emphasized on the point that the maximum density change 
could be about 35 yi • Even though this effect is not too 
significant/ it would not be correct to assume incompressibility 
because a very high pressure in EHD lubricated films is 
encountered* The other significant contributions on this 
effect (of isothermal) of EHD lubricated line contact are 
due to Archard et al [106]/ Dowson and Whitaker [l07j / 

Adams and Hirst [l08]/ Rodkiewicz and Srinivasan [l09]/ 

Cheng [llo]/and Evansand Snidle [ill] etc** 

It was mentioned earlier that the density might change with 
temperature and pressure. This aspect was elaborated by Cheng and 
Stemlicht[ 33-89 j using Newtonian lubricants* They studied the 
very complicated case of thermal EHD line contact and obtained 
pressure/ temperature and film profiles in addition to frictional 
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force. Hirst and Moore [112] also studied the lubrication 
of similar geometry but they assumed the fluid to obey 
the non-Newtonian Maxwell model* They examined this problem 
using different fluid models and elaborated results through 
figures and tables* Apart from this there are several other 
publications regarding compressibility effect in lubrication 
[26/113-114]. 

7. INERTIA EFFECT l 

The classical Reynolds equation is in general derived 
under the assumption that the thin lubricant film is laminar 
and lubricant inertia is negligible* However 'with the current 
trend towards the design procedure for bearing operating with 
high speeds/ and use of lubricants with lower viscosity 
and higher density/ the fluid inertia becomes significant 
for some range of moderately large Reynolds numbers* This 
effect becomes even more significant in case of thick film 
bearing under squeezing motions* 

Smallgpand Vohr [115] studied the effect of inertia in 
hydrodynamic lubrication (of isothermal condition) theoretically 
and experimentally and elaborated the significance of fluid 
inertia* Wada et al [53] observed that in most of the squeeze 
film bearings# film thickness is developed due to the viscosity 
of lubricant* The fluid inertia attains considerable importance 
as the film thickness and squeezing motion increase* Elkouh 
[116] studied inertia effect in laminar flow of power law fluids 
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between two circular plates and concluded that load correction 
factor# defined as the difference of loads (with and without 
inertia ) divided by absolute value of load without inertia# 
increases as n decreases* The non-Newtonian model has also 
been used by Salem et al [69] who incorporated inertia effect 
in a spherical bearing including thermal effects- They 
elaborated pressure and temperature as a function of n* 

Very recently Szeri [98] presented a brief survey on inertia 
effect in lubrication* 

8* IN THIS WORK t 

Classical lubrication mechanics constitutes the solution 
to the standard form of Reynolds equation (i*e* based on 
idealized conditions of isothermal # highly viscous Newtonian 
flow etc*) with various boundary conditions simulating various 
bearing configurations* However much of the modem lubrication 
theory deals with departure from standard form of Reynolds 
equation [6]* It was pointed out earlier that such departure 
may be due to thermal effects# inertia effects# undetermined 
film thickness and various types of non-Newtonian fluid 
behaviour etc*. Most recent information [117] about advances 
in hydrodynamic theory also suggests - the incorporation of the 
following * thermal effects# use of accurate rheological 
models for the lubricant# evaluation of stress. In the materials 
and considerations of unsteady conditions of loading* The work 
in this thesis Is devoted to the study of some of these 
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important factors responsible for the departure • 

The main factor which is responsible for the departure/ 
in heavy loaded/high speed bearings/ is the non-Newronian 
aspect/ which is invariably encountered* During the last 
few years/ power law non-Newtonian model has found extensive 
applications* Keeping in view the motivation of recent 
research and the practical importance of this non-Newtonian 
model/ the power law model is used throughout this thesis- 
In addition/ roller bearing which has advantage over the 
other bearing configurations is chosen for study in this 
work. 

Another factor responsible for the departure is the 
effect of temperature (and pressure) on the lubricant 
consistency m of the power law fluid* Hence for the 
entire work in the thesis/ m is considered to vary with 
temperature (and pressure)* Subsequently tne consistency 
m is assumed to vary with the mean temperature (and 
pressure) [118]'* Use of this mean temperature not only 
helps in reducing the complexity of the problems but also 
yields solutions for temperature in closed form. * 

Other factors considered here are the lubricant compress- 
ibility and the inertia effects* Squeezing motion of the 
rollers with and without inertia has also been incorporated* 
Finally the rollers are assumed to rotate with unequal 
speeds • 
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Table * 1 Representative power law constants 


Material 

m 

n 

54-3°/ Cement rock in water 

25.1 

0.153 

23*3/^ Illinois yellow clay in water 

55.5 

0. 229 

^olystyr ene at 300 °P 

l . 6xic> 6 

0.400 

L*5°/‘ Carboxymatnylcellulose (CMC) in 
atar 

97 

O'. 400 

).7 / CMC in water 

15 

0.500 

5 Polyisobutvl ene (Enjay 

fistannx L-100) in dacalin 

9.4 

0.770 

1-5°/ Separan 

1578 

0.119 

> °/> Hydroxyathyl Celulose in 50*50 
tfatar glycerin mixture 

128 

0.400 

2.5 Weight / polyacrylimide 

Suspension in a 35 I 65 water • 
jlvcerin mixture 

86 

0.545 

Newtonian 

0-75 

1-000 

2th- glycol/H^O/glycerine 
d.th Sp. gravity 1.240 

0.56 

1-150 


Sth» glycol /H^O with Sp. gravity 

..222 


0. 25 


1.180 



CHAPTER I 


THERMAL EFFECTS IN LUBRICATION OF ROLLERS BY A POWER LAW FLUID 
1 . 1 INTRODUCTION l 

The isothermal form of Reynolds equation is commonly 
used to predict the performance of machine components like 
rolling and sliding element bearings* This standard form of 
Reynolds equation is valid for the lightly loaded bearing 
where material properties of the lubricant may be assumed to 
be constant [119-122]. 

However# in case of a heavily loaded bearing system under 
hydrodynamic and EHD conditions/ where pressure is very high/ 
material properties of the lubricant cannot remain constant 
but vary with pressure [l2s]» Hersey and Lodenslager [124] 
reported an investigation of the film thickness between rigid 
gear teeth lubricated by a fluid characterizing a parabolic 
viscosity-pressure relationship. Blok [125] used exponential 
viscosity -pressure relation for analysing gear lubrication 
and concluded that for the fluid following this exponential 
relation/ there was a limiting film thickness at which the 
maximum pressure became infinite. A survey of the earlier 
works in this direction has been presented by Dowson and 
HigginsOn [126]. Recently/ Cameron [127] presented a basic theory 
on lubrication stressing that the lubricant viscosity is an 
exponential function of pressure* Mostofi et al [l28j also 
used the same viscosity-pressure relation to analyse 'the EHD 
lubricated point contact regime at moderate loads and concluded 
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that pressure and film thickness vary with material properties 
(in addition to load and squeeze velocity)* Very recently/ 
an isothermal theory of EHL of rollers having sinusoidal 
roughness/ using the same viscosity-pressure relation has been 
presented by Karami et al [l29j* It was observed that use of 
this relation results in a significant influence of roughness 
amplitude on the distribution of hydrodynamic pressure and 
film thickness at constant load and constant roughness wave 
1 ength« 

It has been pointed out earlier (section 4) that heat 
is generated during the process of lubrication* Imposition 
of squeezing velocity to the bearing surfaces is likely to 
generate considerably more heat inside the system [130-I3l]» 
Because of strong dependence of viscosity on temperature/ it 
is also important to consider the effect of temperature on 
viscosity of the lubricant* Zienkiewicz [132] was perhaps the 
first who initiated the thermal effects in lubrication by 
assuming a temperature dependent viscosity* Dowson and 
Whitaker [133] described a numerical procedure for calculating 
the pressure and temperature distributions and lubricant film 
profiles at heavily loaded/ line contact allowing for the 
consequence of viscous heating in the lubricant* Agrawal 
and Wilson [134] also presented a numerical analysis of the 
flow of a liquid lubricant (characterizing exponential viscosity- 
temperature relationship) under conditions where visco us 
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heating leads to a substantial viscosity variation across 
the xilm. More recently Heshmat and Pihkus [135] used a 
similar model in order to study mixing inlet temperatures 
in hydrodynamic bearings* Giordano and Boudet [136] studied 
the flow of lubricant between two parallel plates having 
unequal speeds/ assuming the lubricant viscosity to vary with 
both pressure as well as temperature* He analysed a THD 
case and on comparing with the isothermal case/ ha found a 
significant difference between these results* The importance 
of thermal effects has also been verified experimentally 
[137-138]. 

In addition/ the tendency towards * limit design* increases 
the need for better tools to calculate bearing characteristics 
accurately. Thus the selection of an appropriate fluid model 
to analyse the operating characteristics in a better way/ is 
also necessary. Although/ Newtonian model with exponential 
viscosity variation with respect to pressure and temperature [ 139 ] 
is the simplest and the most widely used one. But it cannot 
meet the specific requirements of many engineering applications 
where different additives in the lubricating oils are used* 

These oils with long chain polymer additives have non-Newtonian 
behaviour i.e* the viscosity of the lubricant is not constant/ 
but usually exhibits non-linear relation between shear stress 
and the strain rate [l40]. One can therefore/ assume the 
fluid to be non-Newtonian. The power law fluid is one of the 
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most viable and practical non— Newtonian models which has 
been used extensively in recent years [l4l]* 

Despite the tact that the fluid behaviour in lubrication may 
be non-Newtonian and variation of the material properties of 
the lubricant (especially apparent viscosity) is important/ 
there are a very few published information [ 92/ 95--96/142 J 
where these effects have been incorporated*. However# following 
are probably the only publications where consistency (m) variation 
of the extensively used power law model has been accounted 
for ' Salem et al f 69] who analysed spherical bearing 
lubricated by power law fluids (characterizing exponential 
consistency- temperature relationship) / did not consider the 
effect of pressure on m* Sinha and Raj [67] studied roller 
bearing considering consistency variation with regard to 
pressure* Raj [75] further extended his work to include thermal 
effects using the second approach [78]. This approach fails 
to provide any idea about the lubricant temperature field* 

Shukla and Isa [143] considered the effect of temperature and 
pressure on m* However only the squeezing motion of bearings 
was considered. 

The present investigation endeavours to analyse the 
performance of thermal effects on the lubrication of rigid 
cylindrical solid by an incompressible non-thixotropic non- 
Newtonian power law fluid considering squeezing motion and 
cavitation* The modified Reynolds and energy equations allowing 
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exponential .consistency variation along the film thickness* 
are obtained for pressure and temperature respectively and 
are solved numerically* 


1 . 2 MATHEMATICAL ANALYSIS • 

B asic Equations » To obtain an expression for rheodvnamic 
study of lubrication* the Reynolds eguation has to be reworked 
to take account of the abnormal viscosity properties of non- 
Newtonian power law lubricants* The non-linear relation 
between stress and the rate of strain for a power law fluid 
is given [33] by 


T ij ” m 1 /e ij e ji ' n ~ 1 e ij ' ^ 1/2/3 (l ' 1} 

where t_; j is the stress tensor and e . . is the rate of strain* 

-j i J 

m and n arc experimental viscom eerie constants* For a steady 
state lubrication of identical rollers rotating with equal speed* 
as snown in fig* (l.l)* and with the usual assumptions of 
nydrodynamic. lubrication * the flow equations for a power law 
fluid are given [l?9] by 


SjE = 3,1 

ax ay 


where x 


(f rom(l* 


3u 
3 x 

1 ) ) 



takes 


0 

the 


form 


( 1 * 2 ) 


( 1 * 3 ) 


T = ral fy' n ~ X ay 


( 1 * 4 ) 


u and v are the velocity components in the directions of x 
and y respectively and p is the hydrodynamic pressure* 



The most difficult problem in lubricated bearings is 

the determination of th : operating viscosity* This difficulty 

arises from the fact that th.e viscosity of the lubricant 

varies rapidly with pressure and temperature. This variation 

does not lend itself to a simple expression in an analytical 

term* Various theoretical models have been postulated by 

different researchers* In this thesis/ exponential relationship 

for viscosity has been used throughout* Hence one can write 

ftp— 0 (T~T ) /, e ) 

m = m e F o (.JL.5J 

o 

where m Q is the value of m at p = 0 and temperature T = T q 
( ambient) / cc and are the coefficients of pressure and 
temperature respectively* 

The flow is assumed to be symmetric about x-axis* Thus/ 

it is sufficient to consider the region for y ^ 0- The 

boundary conditions for the problem are 

TT . h 

u = U at y s= v? 

. (1*6) 

= 0 at y = 0 

ay 

h is the variable film thickness given [107] by 

h = h Q + x 2 /2R (1.7) 

where h Q is the minimum film thickness and R = r/2/ r being 
the radius of the either cylinder* 

Reynolds Equation • Considering the appropriate sign for 
pressure and velocity gradients [33] for two different 
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regions as indicated in fig* (1*1)/ integration of equation (1.2) 
yields'*" 


u = U " ( m“ l ' /n C(|) (n+ 1 ) - /n -y (n+1) ' /n }/ ~°°£x £ -X. 


u = 


u + (- S- f (§> (n+1) - /n .yfn+lV*,, 

2 


1 

( 1 . 8 ) 


“ n-i-1 


^ X ^ X 2 


(1.9) 


where x = -x^ and x = x 2 are respectively the point of 
maximum pressure and the point of cavitation* 


Integration of equation (1*3) with boundary conditions 


v = 0 


y = o 


(i.io) 


gives [139] 


v = v, at y = 


h/2 

3 f , U dh 

- v h = h J u dY - 5 mr 


( 1 . 11 ) 


where v^ is the resultant of the squeezing velocity V(-| is the 
velocity of the either cylinder) and the normal component of 
rolling velocity U i»e* 

TT TT Jt 

( 1 * 12 ) 


v - + Y u dh 

h “ ± 2 + 2 dx 


Substituting u from equation (1*8) in equation (1*11 ) t the 
Reynolds equation is obtained as 

fe !< 2ifr -) ( 5> (3n ‘ : ' 1)/n <=- ^l 1 /n -a) = V, — s XS - 


m 1 dx 


*1 

(1-13) 


1* m is a function of x-alone because under adiabatic condition 
T is assumed to be a function of x only (see energy equation). 
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A similar expression is found for the other region ire. 


5 * {(—2n — j (^) ( 2n+l)/n 1_ _2) V n + $ = -V/ -x, £ x 9 

dx 1 ^ 2n+l ; k 2' v m 2 ax s 1 2 


(1.14) 


Equations (1.13) and (1.14) are the modified Reynolds 
equations for the two specified regions. 

Pressure Equation • To solve the above Reynolds equations/ 
the following pressure conditions are used [33]* 


Pj =0 at x = -00 


% = p 2 

a Pl dp 2 at x = -x^ h = h 

33T ■ 35T = 0 

dP ? 

p 2 = dx“ = 0 at x = x 2 ' h = h 2 


(1.15) 


Integration of equations (l- 13) and (1*14) using the above 
conditions yields 


= 2m 1 C n (X 2 -X^+2q(X+X 1 ) ) n /(l+X^) 


2^ 2n+l 


(1-16) 


,^2 a -2m 2 c n (X 2 -X 2 -2q(X+X-, ) ) n /(l+X 2 ) 2n+1 


(1.17) 


where 


/4n+2^ n • i n/2R \l/2 _ _ V ( R_^l/2 

— '-U "> V U ' * g ~ TT ''VU'”' 


o o 


U 2h 


( 1 . 10 ) 


( 2Rh o ) 1 ' 7 ^ 1 (2Rh 0 P^ 2 (2Rh Q ) 1 ' /2 
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3p ? 

Further# use of the condition -g-— 
gives 

X 2 r= X x - 2q 

With the non-dimensional scheme 


0 at x = in (1.17) 


(1.19) 


H =|~=1+X 2 / i = 2m c n « / P= a P 
o 

equations (l»16) and (l»l7) can be rewritten as 


( 1 . 20 ) 


== n^d)"/ H 2 "* 1 

dp 2 = / /LJ 2n+l 

d3T = " m 2 ( “ f) / H 


( 1 . 21 ) 

( 1 . 22 ) 


where 


m. = m e 
1 o 


and 


* h-h . = 


P 2 - T 2 


/ rru = m e 
? o 

f = X 2 ~xJ + 2q(X+X x ) . 


( 1 - 23 ) 

(1.24) 


Energy Equation * Heat generation by viscous heating and 
heat transfer from (or to) bearing both serves to create 
large temperature variation in the lubricant film- However/ 
in most cases/ almost all of the work done on the lubricant 
is stored in it as internal energy [139]/ i.e. system may be 
assumed to be adiabatic. This adiabatic assumption also serves 
to reduce the complexity of the energy equation^ 


The steady-state energy equation for an incompressible 
power law fluid under adiabatic condition may be written 
as [ 144 ] 


P 


dE 

dt 


ij 



(1?23) 



32 


where E is •the specific internal energy of the system and P 
is the density of the lubricant* The conduction terms of 
the energy equation have been ignored completely under adiabatic 
condition [143/146]. 


In a generalized thermodynamic approach/ the internal 
energy must be modified. However the basic ideas concerning 
thermodynamics of systems undergoing deformations are still 
in the developmental stage [59]. Besides in non-isothermal 
analysis of power law fluids/ it is customary to use the 
conventional expression for the internal energy [145] • Thus/ 
the energy equation (1.25) for the system under consideration 
and with the assumptions already discussed/ in one dimensional 


case reduces to 


Pcu 

3x ay ay 


(1.26) 


where c is the specific heat of the lubricant which is assumed 
to be constant* Averaging the energy equation (l»26) over 
the film thickness / one can obtain 


Pc 


dT 

Hx 


h/2 


h/2 


/ udy =■= m / 


o 


or 


dT 

dX 


i = ™, <I) n+1 /(SH 2n+1 )- 


I n- 1 (|£) 2 dy 

ay ay 


~oo g x£ -X, 


dT 

^ = rm 2 (-B n+1 /(SH 2n+1 ) 


/ -Xj ^ X ^ x 2 


where 


T x = P(T 1 -T Q )/ = P(T 2 -T Q ) 

S * l+X 2 -2q(X+X 1 ) / y = p/Pca 


(1.27) 

(1.28) 
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Equations (l*27) and (1.28) are the main equations for 
temperature for the two specified regions* 

S olutions • For the first region r - 00 ^ X ^ -X^ / logarithm!' 
differentiation of (1*23) implies 


^ dm^ 
m l 


dp 3 

cDT 


dT. 

dX" 


(1-29) 


Integration of (1-29) after substituting the expressions 
dp-^ dT^ 

for and -^r- from (1*21) and (1.27) results in 


fn^(x) 


/ — oa £ X ^ ~X ± (1.30) 


X 


/ (1 - 


rl, (?) n 


m 


o 


S' H 2n+1 


ax 


Similarly/ for the next region/ the expression corresponding 
to (l-30) is 

1 


m 2 (X) 


-X, 


1 jc /■-, 7r ? dX f ( -j T f ^ " f J dX -j 

-If U_ - l (1- s“ J :)2n+T^ } 


X 


m 


H‘ 


" X 1 


H 


~X 1 ^ X ^ X 2 

Finally the pressure is obtained from equations (l«2l) and 
(1.22) as follows : 


(1.31) 


p (X) 


x m 1 (x)f n d;<: 


.«> H 


I2n+l" 


- / - eo £ X £ -X, 


(1.32) 


-X 


1 m. (X) ? n dX 


p 2 (x) 


H 


2n+l 


X m (X) (-?) n ax 

{ x --yjsfi X 1 S x s x 2 (1 - 33) 
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These expressions for pressure contain an unknown X^ which 
can be determined by using conditions * 


=0 at X = - m and p^ = 0 at X = 

in (1-32) and (1.33)/ and one can obtain an integral 

~ X 1 m, (x)f n dX X 2 rrL(X) (-5) n dX 

T zr - — r £ — Q 

~ 2n+l ~ f H 2n+1 

H -Xj, H 

where X^ and X^ are related from equation (1*19)* 


(1.34) 


equation 


(1.35) 


This integral equation (1*35) can be solved by numerical 
method giving the value of X-^ and thereafter pressure can 
be obtained from equations (1*32) and (1*33)* Similarly 
temperature can be obtained from solving equations (l»27) and 
(1.28) . 

1.3 RESULTS AND DISCUSSION ^ 

Expressions for various bearing characteristics/ (equations 
(1. 21-1. 22) / (1. 27-1. 28)) / are functions of the flow behaviour 

XT TD 1 /2 

index n/ velocity ratio q(= u^"”2h~"^ ' an< ^ "thermal factor 
y ( - p/Pcoc) . For numerical calculations/ the following 
values are used - 

a = 1.16xl0” 9 dyne"* 1 cm 2 > 0 = 6x10*"“ deg c" 1 / 

ft 1 1 m 1 * 

Pc = 2x10 erg cm" degc“ / R = 2.54 cm / U = 400 cm sec / 

h = 10' 4 cm. 
o 

Results are computed for pseudoplastic and Newtonian fluids 
only i.e. for n£ 1« Results for dilatant fluids (n > l) can 
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also be obtained/ provided the initial value o. the consistency 
m corresponding to n > 1 is known* The values of q are not 
straightforward out one may expect q lying between the two 
given rangesl(a) -10 10 / which is valid for steady- 

state lightly loaded condition [l-ts] (b) -Is* q ^ 1/ which 
is valid for unsteady- state lightly loaded condition [l2lj* 

In accordance with our expectation for steady state heavily 
loaded bearing/ q may be assumed to lie in the range * 

— * 1 ^ q St -t-* 1 * 

Numerical Me thod * To evaluate pressure and temperature which 
develop inside the fluid film/ we use equations (1*21-1-22) 
end equations (1*27-1*28) respectively which contain an 
unknown X-^ This X 1 is calculated by solving the integral 
equation (1*35) numerically for various values of n and q/ 
using Simpson rule and bisection method. The point of 
cavitation is then calculated by using the formula given 
in (1*19) and these values are listed in tables (l*l) and (1*2)* 
Finally/ the simultaneous equations for pressure and temperature 
are solved numerically by Runga-Kutta method. The results for 
dimensionless pressure/ temperature and consistency are 
elaborated through figures* 

Pr essure Distribution l The dimensionless pressure p is 
plotted versus X / in fig. (1*2) for various values of n and 
q * It is observed from the figure that for a fixed n / a 
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numerical increase in ~q results in the shift of the 

position of pressure peak towards the centre line of contact/ 

whereas a numerical increase in +q causes the reversal of 

the trend. Physically/ one may interpret that as the 

surfaces approach each other with increasing normal velocity/ 

more pressure is generated* The obtained pressure profile 

as a function of q is quite similar to that of Dowson et al [l2l]* 

The variation of pressure with respect to n is similar 
to that obtainedbv Safar and Shawki [ 63 ]* Pressure 
increases significantly as n increases whereas the position 
of pressure peak does not change significantly* 

Temperature Distribution • The non-dimensional temperature 
T is depicted in fig - . (1*3) versus X for various values of 
n and q and the location of the point of maximum pressure 
is marked by X-^ on each curve • The temperature profile 
shows that there are two temperature dominant regions / 
one is near the inlet and other near the outlet* These 
temperature regions show that heat of convection plays an 
important role there* The existence of the temperature 
dominant regions in the proximity of the inlet and outlet/ 
only conforms to expectations/ since the viscous heating 
which is caused by non-zero velocity gradients exists in 
these regions* The importance of heat of convection in the 
inlet region has been established by Crook [86] and Murch [14-9]* 



The velocity gradient at the point of maximum pressure is zero* 
This causes the temperature to be almost constant in the 
vicinity of that point* The temperature variation with 
respect to n and q is almost similar to that of pressure in 
the region - °° ^ X £ ~X^. In the other region -X^ £ X^ X^/ 
pressure decreases whereas temperature increases* The effect 
of temperature on pressure peak is elaborated in Tables (1*1) 
and Cl* 2)* It is observed from the tables that the effect of 
temperature is to shift the position of pressure peak towards 
the left which is in accordance with that of obtained by 

4 

Cheng and Sternlicht [89]. 

Consistency Distribution • The main aspect of this analysis 
is to study the variation of consistency m under heavily 
loaded conditions* The interesting pattern of this variation 
versus X is clearly demonstrated in fig* (1*4) for various 
values of n and q? In this figure also the point of maximum 
pressure is marked on each curve as For a fixed value 

of q / consistency decreases with n* This decrease is almost arr- 
ested as n reaches the value 0*545* It should be emphasized 
at this stage that for values of n £ 0*545/ a change in q 
fails to produce a significant variation in pressure/ temperature 
and hence consistency • These curves therefore could not be 
plotted. For fixed n and for all q > 0 / the nature of m 
is quite similar to that of pressure (fig* (1*2)) i*e* m 
goes on increasing/decreasing throughout the regions 
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- “> £ X £ -X, / -X, £ X £ X 0 . In this sense the behaviour 
1 ' 1 2 

of in for all q £> 0/ is similar to that of mentioned above 
except near the inlet where temperature dominates pressure* 
Because of the temperature dominance in this small region/ 
m starts decreasing and continues to decrease till the fluid 
reaches the extreme pressure region* The consistency m 
increases with q in the temperature dominant regions/ however 
in the pressure dominant region / m increases as q decreases* 

1*4 CONCLUSION l 

In this chapter/ a generalized Reynolds equation is 
derived for a power law fluid which lubricates two identical 
rollers having the same speed* The effects of pressure and 
temperature (i*e? often encountered in heavily loaded lubricated 
bearings) on the lubricant consistency is taken into account* 

The squeezing motion of the bearing surfaces is also accounted 
for* 

The effect of temperature on consistency is to shift the 
position of pressure peak slightly away from the centre line 
of contact* However/the pressure peak is found always near 
the exit of the film (with and without thermal effects). 

Squeezing motion of the surfaces also has a significant 
effect on pressure peaks/ the points of cavitation (see 
tables (1.1- 1*3) on pressure and on temperature (see fig*(l*2) 
and fig* ( 1 »3 ))* Pressure and temperature increase as squeezing- 
velocity decreases r In addition/ the consistency change/ 
especially in the zone vhere pressure is maximum/ is also quite 
significant. 
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Fig. 1.3 Distribution of T vs. X. 


CHAPTER*" II 

NON-UNIFORM TEMPERATURE IN NON-NEWTONIAN 
FLUID FILM LUBRICATION OF ROLLERS 

2.1 INTRODUCTION ' 

It has been emphasized in the previous chapter that 
extremes of pressures and temperatures are encountered in 
heavily loaded bearings and material properties of the lubricant 
do not remain constant* Attention has been focussed mainly on 
lubricant consistency variation with regard to film pressure 
and film temperature and a number of references in support of 
it have been presented. However / another significant property 
i . e* density also plays an important role in determining 
lubrication characteristics [l29j* The density P of the 
lubricant is also a function of pressure and temperature [150]. 

In fact P increases .with pressure and decreases with temperature 
if the lubricant is liquid. For a gas lubricant P increases 
with temperature* Density P may increase by about 20% at 
exceedingly high pressures [126J* Since the liquid lubricant 
density P increases with pressure and decreases with temperature/ 
the effective density is expected to lie somewhere between 
the maximum increase due to pressure and the maximum decrease 
due to temperature* Therefore the compressibility effect assumes 
a significant importance* 

Dowson and Whi taker [107] examined the isothermal 
lubrication of rigid as wall as elastic cylinders by lubricants 
having pressure-dependant profiles* They observed that the 
+ Accepted by Journal of Tribology/ Trans* A* S*M*E» 
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relevant lubricant properties which were influenced by high 
Pressure were viscosity and density* They made a comparison 
among pressure curves for compressible fluid (with variable 
viscosity) and incompressible fluid (with and without 
variable viscosity) and found appreciable differences among 
-hem. Isothermal SHL problem for the case of two rolling 
cylinders/ at different speeds has bean presented by 
Rodki ewicz and Srinivasan [109] where lubricant has been 
assumed to be compressible throughout the region* Pressure 
and film thickness profiles have been presented for different 
rolling velocities* Evans and Snidla [ill] studied the 
isouhermal EHL of spheres including compressibility effect 
and viscosity variation. They also obtained pressure 
and film thickness profiles* Hamrock and Jacobson [l5l] 
outlined a numerical solution of the isothermal EHL problem 
for line contacts assuming variable viscosity and density. 
Ihis result was found in good agreement with that of Dowson 
and Higgi nson [ 152 ]. Lubrecht et al [153] also studied the 
same problem for isothermal case considering viscosity and 
density variation with pressure but they applied multigrid 
method in order to calculate film thickness 'and pressure 
profiles* Influence of the compressibility of the lubricant 
on the minimum film thickness and on pressure spike has been 
'-Xamined. Experimental result on the compressibility of the 
lubricant has also been reported [26]. 
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In addition/ variation of lubricant density with respect 
to temperature has also been studied [ll4]* One of the most 
important observations made by Cheng [88] has already been 
referenced in which the lubricant density change with respect 
to pressure as wall as temperature was incorporated* 

Mazumdar [ll3] presented a THD solution of oil journal bearing 
where he considered the viscosity and density to be functions 
of temperature and concluded that this consideration leads to 
a small decrease in the load capacity when compared with 
that using constant lubricant properties* Very recently 
Lebeck [154] also considered the effect of temperature on 
lubricant significant properties/ when he examined lubrication 
analysis of parallel slider bearing theoretically and 
experimentally* In his theoretical analysis/ an exponential 
relationship for viscosity and a linear relationship for 
density has been used* Both the results (theoretical and 
experimental) have been found in good agreement* Non-Newtonian 
lubrication of EHD line contacts together with the compressibility 
effect has also been studied by Houpert and Ham rock [105]* 

Gross [155] presented a brief survey (of the early work) on 
thermal effects in lubrication where he concluded that changes 
in density and viscosity along and across the film can result 
in a load capacity for the otherwise non-load supporting 
parallel film. Apart from them there are several informations 
dealing with variations in lubricant properties and it is 
difficult to recall all of them* 



The present analysis focuses on thermal effects in 
lubrication of rigid cylindrical rollers by a compressible 
power law fluid considering cavitation. The properties such 
as density and viscosity of the lubricant are assumed to 
vary with pressure as well as temperature* The modified 
Reynolds and energy equations/ allowing exponential consistency 
and density variation along the fluid film are obtained and 
are solved numerically* The non-dimensional pressure and 
temperature profiles/ load capacity and traction etc* are 
obtained and are compared with those for incompressible case* 


2*2 MATHEMATICAL FORMULATION s 

B asic Equations • Governing equations similar to (1.2) and 
(1*3) including the lubricant compressibility effect can be 
written as follows * 



a£) 

ay ay ay 


2_e^ + = 0 

ax gy “ - 


( 2 * 1 ) 

(3.2) 


where the consistency m and the density P of the lubricant 
are assumed to vary with pressure and temperature as per the 
following relationships [S6/89J 

_ m ap-e (T-T ) 

m=me h o ( ? , 

o 


p ss p 0 P/Ps mE (T-T c ) 

o 


(2.4) 
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where P is the density at the ambient pressure and temperature/ 
p is the bulb modulus of compressibility and e is the 
coefficient of thermal expansion. Assuming that the flow 
of the lubricant through the bearing clearance (shown in 
tig. (l.l)) is symmetric about the x-axis i*e* both rollers 
are identical and moving with the same speed and T is a function 
of x-alone (see section 1.?)/ one may consider fluid flow in 
the region y 2; 0 only. 

Using the boundary conditions given in equation (1.6)/ 
one can solve equation (2*1) and can obtain 


n — tt "i ^-/ n f /h-\ ( n-i-1 ) / tt ( n+1 ) /ii] 

i ~ u ~ { - n ~r t d sr } { V ' -y • } * 


<» g X ^ -X* 


(2.5) 


a 2 = u + ( 5?r ) ( - h~ a v )1/n f <!> (n+1)/ ft_ y (n+i)/aj , 


-X ' g X g X. 


( 2 . 6 ) 


where x = — x* corresponds to the point of maximum pressure/ 
u l /m l an< ^ u 9 ,m 2 correspond the respective quantities in the 


regions - ■» g x g -x' 


md -x* gxg x 2 * 


Making use of the 


conditions (l«10)/ integration of the continuity equation (2.2) 


gives 


h/2 

- p i v h “ p b (p i u i ;d y 

h/2 

/\ 

-Pi v. = rrr- r Pn u, dy - 


i h - ax 


1 i - 


p,u ^ 

1 • dh 

' T“ dx 


( 2.t) 


where v h is the velocity in y-direction at y - h/2 and may 
be written as [67] 


_ U dh 
v h “ ^ dx 


( 2 - 8 ) 
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A further simplification of equation (2«7) on substituting 


Uj from (2.5) yields 


P, ah 


4- f - lUX1 _ 0 (-2 if 1 ™ dP i.\ 1 / n f hs (2n+l) ,/n-j 

cK L _- - — -g-) (-) J = 

-00 £ X £ -X* 

Similarly for the other region/ on can also get 

, n )( _ 1_ ^2,1/n (h^n+li./n n . _ 

HxL 2 + ^ 2 ^ 2n-i-l m 2 cET ; V 9 ; J - U/ 


X* £ X ^ X„ 


(2.9) 


(2.10) 


Integrating equations (2.9) and (2*10)/ using the following 
condi ti ons; 



bp 1 

cix“ 0 #. 

X = ~3C 

P = P* 

/ h = h* 

(2.11) 

One may 

write 



dp x 

dX“ 

= 2m, c (f) n /H 2n+1 / 

I n * 

- 00 ^ X £ -X* 

(2-12) 

dp 2 

dX 

= -2m 0 C (g) n /H 2n+1 / 

^ n ^ 

- 3 ? £ X £ X 2 

(2.13) 

where f 

g 

= H - P*H*/P 1 
- -H + p*H* /P 2 


(2.14) 

(2.15) 


The dimensionless quantities occuring in equations (2.12) 
and (2.13) have already been defined in equations (l»18) and 
( 1 . 20 ) . 
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The boundary and the matching conditions for the pressure 
equations (2*12) and (2*13) together with the cavitation 
boundary conditions are 


P x =0 / X = - oo 


(2.16) 

p x = p 9 ; p 2 = p 2 / x = -x* 


<3p 9 

(2.17) 

P 2 = dX“ = 0 ' X = X 2 


E nergy Equation « Assuming that the fluid is compressible 
and the convection term of the energy equation is dominant 
(see section 1*2)/ then the energy equation reduces to the 
following (see page 706 / ref* [89]) under adiabatic condition ■ 


P 


c u 
P 


dT 

dx 


m 


BY 


(3ji) 
3 y 


where c is the specific heat at 

ir 

is assumed to be a constant* 




constant pressure which 


( 2 . 10 ) 


Integrating the energy equation (2*18) with respect 
to y over the film thickness h and simplifying it/ one may 
obtain 


dT, 

p l c p d)T 


dT, 


P~c 


2 P dX 


2m l c nf> T l H * + ”5 f 3 

2ln 2 c nC _ET 2 f ^ + 

H * H 2n+X 



( a) n r 


-oog X £ -X* 

-x* £ x 2 


(2.19) 


( 2 * 20 ) 


Use of the condition dp 2 /dX =0 at X 
results in 

H 2 = pV/p* 


X 2 in equation (2.13) 


( 2 * 21 ) 
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dT„ 

asr = 0 at x - x, 

# 

where Pg is the density of the lubricant at X = X 2 * H 
Making use of the dimensionless scheme given below/ 
equations (2.12-2.1?) and (2*19-2.20) may be written as 


( 2 . 22 ) 



dp x 

dx“ 


= 5 1 (i) n /H 2n+1 


( 2 . 2 ?) 


dp 


2 


dX 


-H 2 (g) n /H 2n+1 


(2.24) 


= rm 1 (5) n [?+eT 1 ^H%”^ ; Pl' f£( ^l“^o ) ]/(p*H*H 2n+1 ) 

^ _ 

= rm 2 (g) n [5-£T 2 p f H*e~ Ap 2' fe ^ T 2“ T o ) ]/(P^H*H 2n+1 ) 

where 


5 = pT# P*= p */ p Q ' ? - P/(P 0 c p a ) ' e = e/p# A = 1/«P S 


(2.25) 

(2. 26) 
(2.27) 


= H 




Si, = 5 


(2.28) 


5 = -h + pVe-v'vy 


Sij = Si 0 e^2“ ( *2' T o > (2.29) 

Equations (2.23-2*26) are the final foams of Reynolds and 
energy equations giving pressure and temperature distributions* 
These equations have to be solved under the following boundary 
and matching conditions » 

1 • ’‘. 7 ^ 1 / 

k ' J 


1 062 ' 
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Pj. s= 0 / T^ = T q at X = — oa 

P x = P 2 / T x « T 2 at X = -X* (2*30) 

ap 5 

P 2 = 0 = d3T at X = X 2 


The pressure and temperature equations (2.23-2.26) contain 
an unknown X 2 which can be determined as follows [139]* with an 
assumed X 2 / equation (2»23-2»26> can be solved to give 
indepence of p on X. If the assumed X ? is correct* then 
p will vanish at X = X 2 / if it is not/ then another value of 
X 9 is taken and the process is repeated as often as is 
necessary to make p = 0 at X = X 2 to the required order of 
accuracy. Finally / these equations are solved simultaneously 
for pressure and temperature using fourth order Runga-Kutta 
method. 


L oad l The normal force component in the direction of x is 
given [75] by 

h_ X 2 2h X 2 „ 

W = -2 / 2 p dh =s -2 / dx = (“tx^ - ) f X dx dX 

x - » 

The dimensionless load component W v / defined as W x = W v /(2h r> /o 0/ 


X' 


may be obtained as 


W. 


X 


I- ^ rZ2n+l ” 

-oo n 


m X 2 (g) n 

dX - „ — dx] 


-X~ H 


»n+l 


(2.31) 


Similarly/ the y-component of the normal force is 
X 2 

W y = / P dX 

—%DO 

~ X m..X(5) n X 2 m X(g) n 

or w y = - [ t 4m— ax - t - 2 - +r — ax]- 

— oo rl —A in 

where W Y = Wy/( ^ 2Rh o /«) • 
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Traction force 
by 


The traction force Tf may also be defined 



— OO 


which comes out in dimensionless form to be 


T 

X F 


- X rrL ( f ) 

= [ ' 


H 


v,n X„ - ri n 

2 m-CgJ 


2 n 


dX] (2*33i 


where T ? = T^/C-h /2oJ . 


2.3 RESULTS AND DISCUSSION l 


Theoretical aspects of numerically computed results for 
various bearing characteristics are elaborated through figures 
and tables which follow* These characteristics are functions 
of the flow behaviour index n* The index n plays a very 
prominent role in the present context- Results are calculated 
for the following values of n/ i.e* 

0*4 ^ n£ 1.15 

Dor numerical calculation/ the following set of values are 
used • 

-1 -d 

R = 2.54 cm/ U = 300 cm sec / h Q = 5 x 10 ' cm/ 
oc =1.6 x 10 dyne cm and p g = 8.125 x 10 dyne cm 
7 and £/ defined as 7 = |3/(P 0 c and £ = e/{ 3 are dimensionless 
numbers/ values of which are takem as 0.4 and 0.5 respectively. 


In order to study the qualitative behaviour of consistency 


variation of the compressible lubricant/ pressure and temperature 
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must be computed first. This is achieved by solving the 
simultaneous equations (2.23-2.26). These equations contain 
an unknown X 2 which is evaluated by the method described in 
the analysis* Further/ -X* / the point of maximum pressure/ 
is determined from the relation 


X* = 


-Pm aX /ap s + E(I r 


(2.34) 


where p and T are the maximum pressure and the temperature 
r max max 

at X = -x* . The density of the fluid at the exit X = X 2 


is assumed to be equal to that of the ambient'* 


Subsequent to the computation of X“ and X^/ the 
simultaneous eqs * (2*23-2*26) are solved numerically for 
dimensionless pressure p and temperature T by Runga-kutta 
method. This finally yields dimensionless consistency m/ 
using empirical relations given in (2.28) and (2.29). 


The variation of p / T and m with X for various values 
of n are presented in figs* (2.1)/ (2.2)/ ana (2*3) respectively. 
All these curves have one feature in common that variation 
in n does not change the general shape of the profiles* Tne 
compressible and incompressible curves/ denoted by continuous 
and broken lines respectively/ are presented for several 
values of n» Curves in these two cases are hardly distinguishable 
for values of n £ 1- Furthermore/ for n ^ 1/ there is a fixed 
trend between compressible and the corresponding incompressible 
curves / i*e* incompressible curves for p and T lie below 
the compressible curve* Whereas for n ^ 1*15 such a fixed 
trend does not exist'* 
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Pr essure Distribution • The pressure distribution p versus 
X for various values of n has been depicted in fig* (2*1) 
for compressible as well as incompressible fluids* p increases 
continuously in the inlet region and decreases in the outlet* 
Once the pressure peak is attained/ p falls down very rapidly 
to the ambient as the fluid moves towards the outlet. The 
behaviour of p against X for each n is similar to that of 
Dowson et al [l2l]* It is seen from fig. (2*1) that p 
increases with n* This is in conformity with the results 
of Safar and Shawki [ 63 J - One may also observe that for 
values of n s* 1 the p curves for incompressible fluid/ denoted 
by dotted lines in fig* (?*l) / always lie below the correspond- 
ing curves for compressible fluid. A similar trend was 
observed by Dowson ail'd Whitaker [107J for Newtonian fluid* 
However/ for n = 1*15/ this particular trend exists only 
near the outlet followed by reverse trend in the inlet zone* 

Temperature Distribution • The temperature distribution T 
for various values of n is presented in fig* (2*2)* It is 
interesting to note that T increases in the inlet zone* 

This increase is arrested at the point of maximum pressure/ 

X = -X* . T subseguantlv decreases/ rather slowly to the 
outlet. An increase in T in the region -°°5JX <. -X indicates 
that heat of convection is dominant in the inlet region 
[86/ 149] » This increase in T in the inlet region is also in 
accordance with the view point that the dragging action of 
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the faster levers in the high pressure region causes more 
viscous dissipation in the convergence zone and hence results 
in more temperature there* 

T increases as n increases/ this can be seen from 
fLg* (2.2)* An increase in T with n can be physically 
interpreted. An increase in n signifies an enhanced effective 
viscosity* Thus the resistance to the motion would increase/ 
leading to a higher viscous dissipation* 

Incompressible fluid film temperature T is also presented 
in fig* (2*2) (by dotted lines)* It is observed that 
incompressible film temperature is comparatively less than 
the compressible one for all values of n / except when 
n ^ 1*15 * In this case the incompressible temperature T 
in the inlet zone is slightly higher than that for the 
compressible one* The difference is/ however/ not very 
significant/ both for Newtonian as well as non-Newtonian 
fluids* 

Co nsis ten cy Distribu t ion * The main objective of this 
analysis is to study the behaviour of consistency variation 
of a compressible non-Newtonian fluid with respect to pressure 
and temperature. The consistency m varies with pressure and 
temperature according to relations (2*28) and (2*29). 

The general trend of rn versus X for various values of n is 
shown in fig* (2*3)* It is clear from the figure that m 
decreases near the inlet where temperature f (caused by heat 
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of convection) dominates* This decrease continues until 
such time that the pressure p starts dominating. From this 

point and onwards/ m increases rather rapidly/ and becomes 

. 

maximum at X = -X . For heavy loads under consideration/ 

. 

consistency peak near X = -X is nor unexpected. If one 
looks at the definition of in given in ( 2'* 28) and (2*29)/ 
it is obvious that in the high pressure region/ in must be 
high because pressure p developed inside the film through 
hydrodynamic action is more than compensated by T. In the 
region -X ^ X ^ X 2 / m decreases continuously up to the 
cavitation point* It is because/ the dominant factor p 
on which m depends / decreases rather abnormally in this 
region - 

It is also possible to discuss the behaviour of m versus 
n from fig. (2*3). This might throw some light on the 
possible dependence of n on pressure and temperature. It 
is seen that m increases as n increases* For values of n £ 1.0 
there is no significant change in m* It may thus be safe 
to assume that for low values of n* neither m nor n may depend 
upon pressure and temperature- For dilatant fluids no conclusive 
statement can be made* 

The behaviour of consistency in for incompressible fluid 
is also depicted along with compressibility affect* Though 
its qualitative behaviour is quite similar to that of 
compressible one* But there are minor quantitative differences/ 
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especially for values of 1. However* for n = 1.15/ 
compressibility effect seems to be observable* 

Film Rupture • Points of cavitation and maximum pressure 
for compressible as well as incompressible fluid film are given 
in Table (?*l) for various values of n* Mathematically 
these points are related by relation (2*34)* For incompressible 
fluid it reduces to 



If is evident from the relation (2*34) that X 2 ^ ^ 
or H 2 *£ h" [lo7]. it is seen that the difference between X* 
and X 2 for n^ 1 is large* i*e* significant departure from 
incompressible result. One may thus infer that compressibility 
effect is significant for n ^ 1* but for n < 1* this difference 
is not much* i.e* X* &sX 2 . Thus - or n < 1* density change may 
be neglected. This is also in accordance with the physical 
situation that density of a liquid changes only in case of 
high pressure which is true for n ^ 1* 

For compressible fluids* the point of cavitation X 2 
moves away from the canter line of contact (origin 0) as n 
increases* However* in the case of incompressible fluids * 

X 2 moves closer to origin as n increases [33J* 

Load and Traction * Load capacity and the traction force* 
which are the most important bearing characteristics* are presen- 
ted in Tables (2*2) and (2*3)* for various values of n. It is 
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seen from these tables that both the dimensionless load 

components and and the dimensionless traction force 

T_ increase with n / which confirms the enhancement in 
F 

effective viscosity with n j_63j* 

It is observed that; there are no significant differences 
between the values of ^ / Tp for incompressible and 
their corresponding values for compressible fluid when n^ 1 
(pseudoplastic)* For dilatant fluid the compressibility 
effects are more obvious* 

2 * 4 CONCLUSIONS * 

The problem of heavily loaded rollers/ lubricated by a 
compressible power- law fluid is considered* The consistency 
and density are assumed to vary with temperature and pressure- 
It is further assumed that the flow index n does not vary 
with pressure and temperature* Though this is customary [145]/ 
it nevertheless imposes restrictions on the applicaoili uy 
of the analysis* Another important issue that remains 
unresolved is tire variation of the internal energy for power 
law fluids* 

From the present analysis/ it is observed that for low 
values of n there is no significant effect of pressure and 
temperature on the consistency* Even the compressibility 
effects are not significant* It might be thus safe to 
assume that both m as well as n do not vary with temperature 
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and pressure for low values of n* The seme may perhaps 
ba said about the internal energy* However/ for values 
of n > 1/ nothing can be said. 

A careful assessment of these assumptions is highly 
desirabla. Experiment conducted in view of this analysis 
might prompt researchers to investigate further and resolve 
these issues* 
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Table 2.1 


Compressible - Incompressible 


X* 

X 2 

X*' X 2 

m 

o 

n 

0-1712024 

0.8730804 

0.4706875 

0.560 

1.150 

0.3901324 

0.5011706 

0.4725625 

0.750 

1.000 

0.4520124 

0.4826696 

0.4795313 

0086 

0.545 

0-4700024 

0.4825070 

0.4825000 

0128 

0.400 

Table 2.2 * Compressible 

W X 

w Y 

T 

F 

m o 

n 

1.735096 

0.1305279 

1.762750 

0.560 

1.150 

0 * 4525186 

0.2562533 

0.4539417 

0.750 

1.000 

0.1546492 

0.0757574 

0-1552770 

0086 

0.545 

0.03553386 

0.0167051 

0.0355398 

0128 

0.400 


Table 2*3 : Incompressible. 


w x 

w y 

T 

F 

m o 

n 

2.0906140 

1.1104130 

2.117660 

0.560 

1.150 

0.4615181 

0.2424701 

0.4671063 

0.750 

1.000 

0.1553269 

0.0747965 

0.1559035 

0086 

0. 545 

0.0355680 

0.0166580 

0.0356004 

0128 

0.400 




Fig. 2.1 Distribution of p vs. X 
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Fig. 2.3 Distribution of m vs. X. 


CHAPTER III 


THERMAL AND SQUEEZING EFFECTS IN NON-NEWTONIAN 
FLUID FILM LUBRICATION OF ROLLERS 

3-1 INTRODUCTION I 

It is a matter of fact that in a heavily loaded lubricated 
bearing/ high pressure and temperature are experienced- The 
lubricant properties do not remain constant but become pressure 
and temperature dependent- However/ coupling of energy and 
Reynolds equations makes the solution very complicated and 
one has to work out numerically using fast digital computer- 
Moreover/ sometimes convergence of solution is so poor that a 
lot of computer time is consumed and the application of more 
accurate solution becomes less feasible- In view of these 
uncertainties/ one is tempted to make soma extra assumptions 
which may result solutions in closed form or semi-analytical- 
In the first and the second chapters/ energy equation 
was coupled through consistency and density of the lubricant. 
Temperature was assumed to be a function of x-alone. This 
assumption under adiabatic condition may be reasonable for 
qualitative purpose / but for non-adiabatic case it is not 
accurate because the temperature gradient across the fluid 
film is not small [l3l]* If one relaxes this assumption / 
solutions of the coupled Reynolds and energy equations become 
complicated- Analytical solution is not possible unless one 

^Published in Wear volume 119(1987) 175-190- 
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makes soma other reasonable approximations* 

In this regard/ Radzimoviskey [ 2 ] mentioned that 
due to mathematical difficulties which arise when analysing 
such a bearing/ the influence of temperature change and 
pressure change along (and across) the fluid film is generally 
ignored and viscosity corresponding to the mean film temperature 
and pressure is assumed in the computation* Cheng and 
Stornlicht [89] solved THD lubrication of elastic rollers 
where viscosity variation with film pressure as well as the 
mean film temperature was mentioned* Recently Pascovici [156] 
applied the exponential viscosity variation with the mean 
temperature of the lubricant and compared the theoretical 
results with those obtained experimentally [157]. Both 
the results were found in good agreement* More recently 
Malik et al [158] presented an isothermal analysis of % 
h vd ro dynami c journal bearings using non-Newtonian lubricants 
by Viscosity averaging across the filirw The performance 
characteristic results obtained in this analysis for a plane 
journal bearing with lubricant following cubic shear strain 
rate law/ compare very well with the experimental and 
theoretical results available in the published literature* 

Very recently Hashimoto and Wada [ll8j examined turbulent 
lubrication problems of a sector-shaped til ting-pad thrust 
bearings with thermal and elastic deformations of a pad 
using the exponential relationship for viscosity-mean 
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temperature* Lubrication and energy equations have been 
solved simultaneously giving film pressure and temperature 
distributions* thermal and elastic deformations and inclinations 
of a pad. 

Another reasonable approximation to avoid mathematical 
complication is to neglect the heat-convection of the enercrv 
equation in comparison to the heat of conduction across the 
film . Crook [86] pointed out that except for a region 
far removed from the pressure zone* convection plays little 
part in dissipation of the frictional heat and the temperatures 
in the pure rolling are effectively those needed to dissipate 
the heat by conduction across the film. Johnson and 
Greenwood [92]/ and Heyes and Montroso [159] studied the 
influence of temperature on viscosity in the thermal analysis 
of non-Newtonian lubricant under EHD line contacts . They 
too had neglected the convection term in the energy equation- 

Thus the present analysis deals with the problem of 
two inelastic heavily loaded rollers* lubricated by non- 
Newtonian power law fluids under rolling and squeezing 
motions* considering cavitation where both the approximations 
mentioned above have been made i«e* the consistency of the 
lubricant is assumed to vary exponentially with film pressure 
and the mean temperature and the convection tern of the 
energy equation is dropped . Subsequently the Reynolds 
equation and the energy equation are solved simultaneously. 

This results into a closed form solution for the temperature 


field. 
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3*2 MATHEMATICAL FORMULATION 


Bas ic Equat i ons * Referring to the previous section / two 
kinds of approximations have been made » (a) viscosity of the 
lubricant varies with the pressure and the mean film temperature/ 
and (b) convection term of the energy equation may be dropped 
in comparison to the heat of conduction across the film. 

Keeping in view these approximations/ energy equation for an 
incompressible power law fluid (with the usual assumptions) 
may be written as [ 32/ 15 6] 


v 9 T 
9 Y 


+ ml 


IS I n> 

9 Y 


9 Y 


= 0 


(3-1) 


where k is the lubricant conductivity and the consistency m 
characterizes the following exponential relationship * 


«p-(3 (T -T ) 
m =5 m ~e mu 

o 


Here T is the mean film temperature defined as 

m 


T = i / 
m n 




>15 


T dv 


(3-2) 


(3-3) 


The other quantities occuring in the above equations together 
with the other governing equations have already been defined 
in section (1-2). 

The boundary conditions for the energy equation ( 3 - 1) are 

— = 0 at y = O (3.4) 

3 Y 

T = T h at y = + ^ 


(3.5) 
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where T^ is the surface temperature which may be assumed 
to be a constant (isothermal boundary). 

The first boundary condition ( 3 « 4 ) may be called 
the symmetric condition (about x-axis) for temperature- 
Using the symmetry of the flow about x-axis and the 
condition (3*4)/ the relation given in (3*3) can be rewritten 
as 

_ h /2 

T = 4 / T dy ( 3 . 6 ) 

m h 0 2 

S olutions • It is obvious from the relations (3*2) and 
(3*3) that m is a function of x only -Therefore the Reynolds 
equations (for the two different regions) remain the same 
as those given in section (1*2). Hence the energy equation only 
will be reconsidered. 


Since the fluid flow is symmetric about x-axis/ the 

velocity gradient can be evaluated by integrating equation 

0 X 

( 1 . 2 ) once/ using the symmetric condition given in ( 1 * 6 )/ 
one can obtain 


aji = ( 2 

3 y m 



- 00 ^ x ^ 


-*i 


3_u 

3y 



d P 


2 -* 1 /n 



x £ x 2 


(3*7) 

(3.8) 


Integration (twice) of the equation (3*1) / using 
the conditions given in (3*4) and (3*5) and relations (3 *7-3* 8)/ 


yiel ds 



m^n 

T h + k(2n+lM3n+l) 


/ 1 ^1 \ ( n-f 1 ) /n r/h%(3n+l)/n (3n-rl)/n] 

''m, 3x“ ( 2 ' ~ Y 

1 (3-9) 
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T - T + -r-H^L- r (- L_ ^2) (n+1) ^f ( h)(3n+l)/n (3n+l)/n, 

l 2 i h + kCTn+niSn+lT m_ ST J V Y 

2 ( 3 . 10 ) 

where m^ / T^ and m 2 / T 2 correspond to the respective 
quantities in the regions - 00 £ x £ -x^ and -x-^ £ x £ x 2 » 

The mean temperature in these two regions can be obtained 
by substituting the values of temperatures T^ and from 
the expressions (3.9) and (3*10) into (3*6). Thus 

m - . J!l n2 __ ( 1_ ^i, (n+1) - /n (h)(3n+l)/n 

T m, T h + kX2n+l) (4n+I) ^m, dx ; 2 


<£ x £ - x , 


(3-11) 


T = T, + 
m 2 h 


! d-P 2 (n + l)/n h On+U/n, 

IT ( “ Sg cET } K i 


-x^ £ x £ x 2 


( 3 . 12 ) 


dp x dp 2 

Substituting the expressions for c[x~ an d ^x” ^ TOrri 
and (1.17) into ( 3 . 11) and (3*12)/ one can write 

V = T h + ( "Sr )m i k i c n Fn+1/H2n - -~ s xS - x i (J - 13) 


V T h + - -X lS XS X 2 

Uh h_ *1 /o 


( 3 * 14) 


Uh n , /9 , 

where \ ( 2 |)' * * = ? (see section (1.2)) and 


G = -P. 


Again using the same non-dimensional scheme as given in 
section (1.2) together with 


T /k, etc. 

ITT 1 


(3.15) 


K = pk. 
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equations (3»13) and (3*14) can be rewritten as 


\ = 5 h + ‘jH+f S 1 J' n ' , ' 1 ./H 2n ' -» s X S -X, 


(3.16) 


5 h + ( 4S72 ) S 2 G n+1 /H 2n - -5^ S X S X 2 


(3.17) 


, - - p - k(t -T ) 

where m = m o . . m o etc* 


(3.18) 


Finally/ differentiation of the above equations (3*16) and 
(3*17) with respect to // X // / gives 2 


? (- s — )m. f m, — *— 

m 1 4n+l 1 L 1 H 2n 


- + 2 { H(n+1) (X+q)~2nXF}] F 1 


H[ + (-2-) IE, F n+1 ] 

u 4n-M 1 J 


-<» £ x £ -X, 


(3*19) 


pU TJL 

fH 2 ^2xT + 2 { H (n+1) (X+q) +2nXG } G 1 


4H 2n + (^r> ^ 2 0 n+1 ] 


-X, £ X ^ X, 


(3*20) 


The equations (3*19) and (3.20) are ordinary differential 
equations having initial conditions 


T = T at X = - co 

o 

T = T m at X a -X, 

nij 1 


(3-21) 


respectively/ which are coupled with the Reynolds equations 
(1*21) and (l.22) viz* 


F n /H 2n+1 


(3.22) 


ux, 

2. = 0 is assumed for isothermal boundaries* 



7 ? 


= -m 2 G n ,/H 2n+1 (3'. 23) 

The pressure conditions for the last two equations are the 
same as those given in section (1*2) (see condition (l*34))» 
Equations (3*19)/ (3* 20)/ (3*22) and (3*23) are the required 

equations for temperature and pressure distributions which 
have to be solved simultaneously using the same technique 
given in the last of the section (2*2)* 


Load and Traction * The non-dimensional load and the 
traction force T^, / as defined earlier/ may be obtained 
as follows * 


_X 1 Xm,F n 

r 1 


-oo H 


2n+l 


dX - / 


^2 Gn 

TT 2n+l 


(3*24) 


~ X 1 m, F n X 2 fnG n 

/ ~~h— dx - r ”2F“ « 

-oo -X x HT n 


(3*25) 


where m 2 are calculated using the relation (3*18). 


3*3 RESULTS AND DISCUSSION * 

The present analysis is devoted to the study of 
consistency variation in a steady-state thermal lubrication 
of heavily loaded cylindrical roller bearings* The semi- 
analytical solutions for pressure/ temperature and the 
frictional force generated inside the fluid film are obtained* 
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Parametric Values • All -che bearing characteristics/ which 

specify the lubrication behaviour of the system/ depend 

upon the various parametersin (flow behaviour index)/ 

q (squeezing ratio defined as r? ^ ~~t ) and X (thermal 

pUh^' 72 ^ 2h o 

factor defined as — — ) « n is assumed to vary from 0.4 

4ka ^ 2R 

to 1*15. The values of q for the system under consideration 
(steady-state heavily loaded) are chosen in the range 
-*1 £ q£ -Hi . The thermal factor K is introduced to study 
the behaviour of lubricant consistency with temperature 
and the value of xvhich is chosen to be 0.5* A particular 
value of K (=o) corresponds to no effect of temperature on 
consistency [67] and a = o/ K = o correspond to the case 
of constant consistency throughout the region [66]* For 
numerical calculation/ the following set of values are 
used * 

-1 , -4 

R = 2*54 cm / U - 400 cm sec • ’ n Q = 6x10 cm / 

-9 -1 2 

a = 1.6x10 dyne' cm * 


Numerical Calculation l In order to study the qualitative 
behaviour of the consistency in which is a function of pressure 
p and the mean temperature T / the simultaneous equations 
(3*19/3* 2*2) and (3*20/3*23) have to be solved first/ to 
compute p and T m * These equations contain an unknown X 1 
which is evaluated by an iterative process (see section 3*2) 
using conditions p-^ = o at X = -“ and = o at X = = X^-2q* 
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Having found and hence X ? / tha resulting simultaneous 
equations (3*19/3*22; and ( 3* 20/3 • 23) are solved numerically 
for p and T m using the Suler-modified method* Subsequently 
m is calculated from the empirical relation (3*18). 

The results thus obtained are elaborated through 
the figures and table which follow* All these curves 
have one feature in common that variations in n and q do 
not change the general shape of the profiles* Moreover/ 
for all n< 1/ these profiles hardly vary/ thus only a few 
values of n and q are considered for these profiles* 

Pressure Profiles • Curves for pressure p versus X for 
various values of n and q are presented in fig* (3*1) and 
a few related results are given in the table (3*1)* It is clear 
from the figure that for fixed n and q / p starts rising 
from X = - oo and continues to do so till it reaches the 
point X = — X-^ • Beyond this/ it starts decreasing with a very 
steep slope and reduces to the ambient pressure p = o at 
the point of cavitation X = 

For a fixed q/ p increases significantly with n / 
especially for n £ 1/ which is in accordance with Safar 
and S'nawki [63J * . From the table one can observe that for a 
fixed q/ the points of maximum pressure and the points of 
cavitation both tend to go nearer to the centre line of 
contact/ (the origin 0) as n increases [67]. For a fixed 
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n# p increases remarkably as q decreases and the cavitation 
points move slowly towards the centre line of contact as 
q increases [66] ♦ This significant change in pressure with 
respect to q accounts for the observation that as the 
surfaces approach each other with increasing velocity/ 
comparatively more pressure is generated* The pressure 
curves presented in this analysis for various values of q 
are similar to those obtained by Dowson et al [l2l] for a 
Newtonian fluid* 

Temperature Profiles » The influence of n and q on T m 
has been illustrated in fig* (3* 2)* It is clear from the 
figure that for a fixed value of n (= l) and q / T m decreases 
throughout the inlet region up to the point X = -X^/ then 
it increases upto the origin 0 /and thereafter decreases 
up to the point X = X^. However / increase or decrease in 
T m beyond the point X= -X^ is hardly discernible and may 
be treated as constant- The same trend is true for n <1 
also though it has not been shown in the figure* The 
Newtonian curve (n=l) for T m / shown in fig* (3*2 )/ is almost 
similar to that obtained by Crook [86] / T m follows a similar 
trend for n > 1 and q> o* However/ for n > 1 and q £ o/ 

T m increases initially (though only marginally) then decreases 
and finally follows the Newtonian behaviour. This trend 
is almost similar to the experimental result obtained by 
Pascovici [l57j • 



Letting n 1) fixed and varying q in the specified 
interval (-»1# +*1J# T decreases as q increases near the 
inlet while in the proximity of the surfaces T m increases 
with q« However for a fixed n > 1/ T decreases as q 
increases throughout the region except in the proximity 
of the surface (i-e* in the range/ -7 S X £ X 0 ) where T m 
increases with q» The behaviour of T with n can be inter- 
preted from fig- (3*2) * It is observed that T m for Newtonian 
(n=l) and non-Newtonian are not significantly different 
in the inlet region/ however/ as one proceeds towards the 
point X = -X-^ / the behaviour of T m is markedly different 
for both dilatant (n > l) as well as Newtonian fluids* 


The mean temperature T for K = 0 is also calculated 

l 

for various values of n and q* Since the general behaviour 

of T with respect to n and q remains the same as that of 
m i 

T m / so only selected values of n and q have been considered 

for T - From fig* (?*3) / it is clear that for n £ 1/ 
x 

T is less than T m except in the n'hood of the point 
i 

X = — X*i » For n > 1/ values of T near the inlet .are. larger 
1 ■ m l 

(though only marginally) than the corresponding values of 


T_ 


m 


While in the remaining portion/ T is less than T 


m. 

x 


m 


Over all T is quantitatively less than T for all n and q- 
m i 

One may /thus# conclude that the assumption of temperature 
independent consistency leads to under*- estimation of temperature* 
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Consistency Profiles • The curves for consistency in are 
presented in fig. (3*4) for various values of n and q. 

For fixed values of n(^ lj and q / qualitative behaviour of 
m is identical to that of pressure* The same trend is 
followed for n > 1 and q > 0. However/ for n >1 and 
q ^ 0 / m decreases initially as a consequence of increased 
temperature T^ near the inlet. Thereafter/ in increases 
and follows the pressure (p) trend. This is an immediate 
consequence of the profiles for p / T m (given in figures 
(3.1) and (3*2)) and the empirical relation (3*18). For a 
fixed n / m increases with q near the inlet while in 
the proximity of the surface m decreases as q increases* 
However/ for a fixed q/ m increases with n. It is interesting 
to note that significant changes in p and T m which appear 
in the neighbourhood of the point X = -X^ / cause a correspond- 
ing change in consistency. It is also noted from the table 
and figure that the effect of temperature on consistency 
is (i) to shift the position of pressure peak and the 
cavitation points slightly towards the centre line of contact 
and (ii) to make the slope of m very sharp. The change in 
consistency for dilatant and Newtonian fluids/ especially 
in the neighbourhood of the point X = -X^ / is very remarkable 
which justifies the inclusion of consistency variation in 
the treatment of the lubrication of a heavily loaded system. 
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The consistency in. (for K = 0) is also computed for 
various values of n and q. It was seen that the behaviour 
of in relation to in was similar for all n and q/ hence 
only one value of q has been chosen for the graph of rru 
(fig*(3*5))* For all n(S l) and q* rrn is quantitatively less 
than the corresponding values of m throughout the region* 

The same trend is followed for n > 1 except near the inlet 
where rru S m (though only marginally)* 

Load and Traction * The non-dimensional load component 
W for the system under consideration is presented in the 
table (3*1) for different values of n and q* It is observed 
from the table that the load component W increases with 
n for a fixed value of q [63/66]* However/ it decreases 
as q increases for a fixed n, which is expected/ since p 
increases as q decreases* 

The load component W (for K = 0) is also presented 

Y i 

in the same table for all the considered values of n and q* The 

qualitative behaviour of W for all n and q is also the 

Y i 

same as that of • However/ W is quantitatively less than 

x i 

W^» This follows from the fact that value of consistency 
m is higher than that of in. / so is pressure and thus 


the load : 
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The non-dimension ali zed traction Tp is also 
presented along with load component in the same table* 

Its behaviour with respect to n and q is the same as 
that of load component for both the cases K = 0 and 
K + O* 

f 

? * 4 CONCLUSION l 

The problem of pure rolling of two inelastic identical 
cylinders lubricated with an .incompressible power lav; fluid 
is considered including squeezing motion and cavitation* 

The heat of convection which is small as compared to that 
of conduction is neglected* The lubricant consistency is 
assumed to vary exponentially w T ith pressure and the mean 
temperature, Reynolds and energy equations for pressure 
and temperature are derived and are solved numerically 
for various bearing characteristics with and without 
consistency variation with temperature* 

The effect of temperature on the lubricant consistency 
is to shift the position of pressure peak slightly towards 
the centre line of contact* However/ the cavitation goes 
away from that contact* The flow index n and the squeezing 
parameter q are found to have pronounced effects on various 
lubrication characteristics such as pressure/ the mean 
temperature/ consistency / load and traction (see figures and 
table) * The consistency variation/ particularly in the 
pressure peak region/ for the Newtonian and the dilatant fluids 
is quite significant* 
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Table 3*1 J Comparison of bearing characteristics 




Thermal : 

factor K / 0 



n/m 0 

q 

X 1 

X-2 

w y 

T 

F 

1.15/0.56 

-0.09 

-0-05 

0.00 

+0.05 

+0.09 

0-389453 

0.415625 

0-449420 

0.484844 

0.514453 

0.569453 

0.515622 

0.384844 

0. 334453 

2.49206 
2-30833 
2.10 661 
1-92642 
1.79434 

1.79640 

1.55546 

1.31386 

1.11804 

0.93627 

1.0/0.75 

-0.09 

-0-05 

0.00 

+0.05 

+0.09 

0.413750 
0.436875 
0.468125 
0. ^02187 
0. 529375 

0.593750 

0.536875 

0.40 2187 
0.349375 

0.363025 
0.348435 
0. 330036 
0.311594 
0.297275 

0.215015 

0.197837 

0.177845 

0.159402 

0.146149 

0.545/86 

-0.09 
-0.05 
0.00 
+0 .05 
-HD. 09 

0.423750 

0.447187 

0.478750 

0.511719 

0.540312 

0. 603750 
0.547187 

0.411719 

0.360312 

0.120024 

0.117723 

0.114665 

0.111471 

0.108768 

0.063030 

0.059685 

0.055638 

0.051769 

0.048802 

0.4/12.6 

-0.09 
—0 *05 
0.00 
+0 »05 
+0-09 

0.426250 

0.451250 

0.485000 

0.510750 

0.545000 

0.606250 

0.551250 

0.418750 

0.365000 

0.027 369 
0.027009 
0.026503 
0.025970 
0.025534 

0.013770 

0.013137 

0.012374 

0.011644 

0.011091 


Continued 


Thermal 


n/m 0 q 


-0.09 0.416797 

-0.05 0.439922 

1.15/0.56 0.00 0-470703 

+0.05 0.503906 

+0.09 0. 53l 7 19 


-0.09 0.418750 

-0.05 0.442500 

1.0/0.75 0.00 0.473750 

+0.05 0.504370 

+0.09 0.532500 


-0.09 0-424687 

-0.05 0.448594 

0-54/86 0.00 0.480625 

+0.05 0-513750 

+0-09 0.541562 


-0.09 0.429150 

-0.05 0.453 875 

0.4/126 0.00 0.488125 

+0.05 0.519500 

+0.09 0.54 6500 
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factor K = 

0 


X 2 

W Y. 

X 

T 

F . 

X 

0.596797 

2.22211 

1.45908 

0*539922 

2.05843 

1.27704 


1.877 81 

1.08956 

0.403906 

1.71653 

0.93476 

0.351719 

1.59946 

0.82987 

0.598750 

0.352244 

0.205313 

0.542500 

0.337798 

0.188885 


0.319833 

0.169906 

0.404370 

0.302528 

0.152991 

0.352500 

0.288468 

0.140203 

0.604687 

0.117500 

0.0612635 

0.548524 

0.115256 

0.0580430 

— 

0.112295 

0.0541050 

0.413750 

0.109188 

0.0504130 

0.361562 

0.106607 

0.0475800 

0. 609150 

0.027196 

0.013661 

0.551875 

0.026860 

0.013044 

— 

0.026381 

0.012000 

0-419500 

0.025846 

0.011577 

0.366500 

0.025426 

0.011037 
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Fig.3.2 The mean temperature Tm vs.X for K#0. 



Fig. 3.3 Comparison between T m and T mi 


-4 —3 -2 -1 0(0.0) 1 X 

Fig. 3.1^ The consistency distribution m vs. X for K # 0 . 



Fig. 3.5 Comparison between m and m,*. 


CHAPTER IV 


THERMAL AND INERTIA EFFECTS IN HYDRODYNAMIC 

LUBRICATION OF ROLLERS BY A POWER LAW FLUID 

4*1 • INTRODUCTION t 

In the preceding chapters/ various aspects of thermal 
effects in lubrication of roller bearings have been studied* 
Lubricant inertia was assumed to be vary small in comparison 
to the viscous force* However/ to meet the present day 
requirements in science and technology/ bearings are subjected 
to a variety of situations such as high operating speeds/ 
the use of lubricants with low viscosity and high density etc* • 
Under tnese circumstances fluid inertia is likely to be 
significant [161/162]* This effect is important even in 
case of laminar flow [163— 165 ] where an increase in load capacity 
is observed [166]. Inertia forces attain greater significance 
if the bearing surfaces are also subjected to squeezing 
motion [167]. An increase in hydrodynamic film thickness 
may also sometimes compel one to consider this effect* 

Milne [168] conducted an order of magnitude analysis of the 
terms in the Navier-stokes equations* It was observed that 
the influence of the lubricant inertia is most prominent 
in the region of large clearence and for rapid changes in 
velocity* Ting and Mayer [169] investigated inertia and 
thermal effects simultaneously in a parallel surface s<-ep 
thrust bearing- They found a correlation between their theory 
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and the experiment performed by Coombs and Dowson [l7o]* 

Smalley and Vo'nr [115] studied lubricant inertia for a very 
arbitrary geometry and demonstrated convincingly its significance 
in lubrication- Dowson et al [l7l] conducted a wide range of 
experiments using a cylinder on a plane* It was observed 
that/ for a moderately large Reynolds number/ the experimental 
data deviated substantially (due to inertia effect) from 
the Coyne-Elrod rupture conditions* Very recently You and 
Lu [l72] presented a cylinder-plane bearing and a journal 
bearing analysis theoretically/and showed that the 
lubricant inertia had a profound effect even at moderate 
value of the Reynolds number* 

In recent years researchers have directed their efforts 
to analyse inertia effects in fluid film luorication using 
non-Newtonian fluids* Elkouh [116] analysed the effects of 
inertia forces in non-Newtonian squeeze films using power law 
fluids* It has been shown that inertia correction in load 
capacity is more significant for pseudoplastic fluids (n < 1) • 

In another paper Elkouh [173] used the same non-Newtonian 
model to analyse this effect in squeeze film between two plane 
annuli and presented various other bearing characteristics* 

The same non-Newtonian model was also used by Salem et al [69j 
who presented a spherical rolling bearing analysis including 
inertia and thermal effects simultaneously* They illustrated 
frictional power as a function of n in addition to pressure 
and temperature distributions* Tichy [174] investigated 
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circular plates -squeeze film using non-Newtonian viscoelastic 
lubricants and elaborated the effect due to inertia on various 
bearing characteristics* Hashimoto and Wada [53] also 
investigated theoretically the effects of fluid inertia 
forces in parallel circular squeeze film bearings lubricated 
with pseudoplastic fluids and observed that these effects 
become considerably significant with increase in squeezing 
velocity as well as the film thickness. 

In view of the present specific requirements where the 
bearings operate with high speedsand under extreme conditions 
of pressure and temperature/ the consideration of inertia 
effect seems to be important. In the work presented here/ 
inertia and thermal effects in a squeeze film roller bearing 
lubricated with a power law fluid are examined theoretically* 
The fluid consistency m is assumed to vary with pressure 
and the mean film temperature* The modified Reynolds and 
energy equations obtained here are solved simultaneously* 

4*2 MATHEMATICAL FORMULATION l 

Basic Equations * Making use of the constitutive relation 
for a power law fluid (given in (l*4) )' applying the usual 
assumptions of hydrodynamic lubrication for the geometry shown 
in fig* (4*1) and retaining the fluid inertia term/ one obtains 
the following equation of motion [116] * 
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P( u + v + ^ = n- ' i ‘ ^Ji) (4*1) 

ax ay dx ay ay ay ^ 

The first term in the left hand side of this equation is 

the fluid convective inertia which is usually neglected in 

hydrodynamic lubrication. 

The equation of continuity is 


m = o 

ax ay 


(4.2) 


The temperature -within the lubricant film resulting from 
viscous dissipation is defined by the energy equation [143] 


Pc(u fx + v ay^ 


,a_U|n-l ( iU)2 


(4*3) 


Boundary conditions for these governing equations are 
considered to be the same as given earlier by equations (1.6)/ 
(l»10)/(3»4) and ( 3 • F ) * The relation connecting m/p and T 
is also taken to be the same viz* 


a p -p(T -T ) 
m=m e c r m o 
o 


with 


T = r 
m h 


h/2 

/ T dv 


(4.4) 


(4.5^ 


Making use of the folio wine non-dimensional scheme • 

Y = 2y/h / u = u AJ / v = ~ (— ) X ' /2 / 5 = pT 

o 

_ Ch h -l fry 0 

X = f>u4 , P r = jgSjj (,2_)1. , E t = uVc 


(4*6) 


Equations (4 *1-4.5) reduce to 

X(5|| + v||) + §§ = A^m i-l-l 11 " 1 |~) 


(4.7) 


+ IZ _ 0 
dX + 3Y ~ U 


(4.8) 
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if + V ||) = -I- m P„ Ej_A iyi n-1 (-|^) 2 


9Y 


3Y 


r t; 


3Y 


3Y 


p-(T -t ) 
m = m m o 

o 


1 


H 


T = i / T dY 
m Ho 


(4*9) 

(4.10) 

(4-11) 


where A = ( n/( 2n+l) ) n • The other quantities occuring in 
equations (4.7-4*11) have already been defined in chapter I* 

The presence of inertia and convection terms in equations 
(4*7) and (4*9) respectively makes them complicated • Therefore 
theiy” average values are considered. Consequently equation 
(4*7) may be written as 


L / (uiu -2i )dY dg_ A § — ,aii| n -l 2H) 

H Q d X + V 3Y ;ax + dX ~ A 3 Y Vm 9 Y dY J 


( 4 * 12 ) 


- 3 T 3 

The contribution of the second term v is small as compared 


**** T 

to u f-r? and therefore it is neglected from (4*9). Further if 

dT 
n 

ax ~ dx 


one assumes ^7 & [89]* Equation (4*9) now becomes 


dT_ 

K P_ U m -, v 
r m dX 


m 9 2 T , - p i SLH 1 ri— i Aun 2 

“2 + m P r E t A 'aY 1 ( §T 

dx 


(4.13) 


where the velocity u m is taken as 


Letting 


= ^ { y=o + u1 y _hV 2 


f = s f o (a B + ^ l§ )dY + H 


ai 

5 = ^ p r S m ax 


m 


( 4 . 14 ) 

(4.15) 
(4* 16) 


j~» v i*e* u yy >> v -9^ / see ref [89]/page 702. 
ax ay 9* a y 
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It may be noted that F and D are functions of x-alone« 
Solving equation (4*12) for u / using boundary conditions 
given below ' 


§Ji - o 


v = o 


u -- 1 


at Y = 0 


(4-17) 


o , dH 
v = -2q + ^ 


at Y = H 


One obtains 


u = 1 + (-—■) (—£—) 1 / n ( Y ^ n+ l ^ / n _ H ( n+ l ) / n ), _ oo ^ x ^ -X, 
n+JL * ~ i 

( 4 » IB) 


Am, 


The volume flux Q / defined by 

H 


Q = / u dY 

o 


(4-19) 


is obtained as 


Q = H- ( ~- - r ) ( — — ) 1>/n ^n+l) /n ( 4 . 20 ) 


'2 n+l 


Am, 


Integration of continuity equation (4*3)/ using the boundary 
conditions (4*17) gives 

H 


3X 


/ u dY = -2a 


( 4 . 21 ) 


or Q = -2qX + d 

where d is the constant of integration. 

Use of equations (4*18)/ ( 4 . 20 ) and (4*21) in equation ( 4 . 12) 
yields 

dp, 

(4.22) 


I n (X) + d 3T = J f (X) 
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wnere 


- -n /t r2n+l 


f = H+2qX-d * J^(x) ^ m^f/H 


I (x) = £ [2X£ -1 + 


2B f B f 


n 


B f 


n 


H 


H 




with 


B n = 2 ( 2n+l) /( 3n+2) 
Now |“ = 0 at X = -x 1 3* 

f = X 2 - X 2 + 2q(X+X 1 ) 


(4* 23) 


(4*24) 


where the unknown X^ is positive- 

Similarly equation (4*13) can be solved for T- Using the 
boundary conditions 

^ = 0 at Y = 0 

(4*25) 


3Y 


T = T h at Y = H 


the solution is obtained as 

- y“ -F n+ 1 

T i = 5 h + §(^ 2 -h 2 ) + vferiTTwaT (H (3n+1) - /n -Y <3n+1) / n ) 

hi — - 

(4.26) 


n 


where r = 

Averaging (4*26) using definition given in (4*11) / the 
expression for the mean temperature is obtained as follows 


T, 


m. 


-h 


o 

H rr , „ - j=n+l r ,2n 
—— D + r, i, x /H 
3 11 ' 


(4-27) 


where 7 ± = * 
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Substituting the value of D from equation (4*16) into (4*27)/ 
one gets# 


2 d V 

X P r 3~ dx~" = T h 


m. 


v - -n+1 ,„2n 
r, m. f /H 
1 1 


(4.28) 


The differential equations (4*22) and (4*28) are the final 
forms of Reynolds and energy equations giving pressure and 
the mean temperature distributions for the first region 
£ X £ -X^ # where f is given by (4*24) and u (using 

(4.14)) is aiven below* 


n - 1 _ (££±£) 

m v n+l' 2H 


(4.29) 


Similarly# for the other region# -X^ ^ X ^ 0# equations 
corresponding to ( 4 . 22 )# (4*26) and (4*28) giving pressure# 


temperature and the mean temperature distributions are 


dp, 

J n (x) + cDT 


J_ f (X) 


( 4 . 30 ) 
(2n+l) (n+l) 


T 2 = T h+ §(Y 2 -H 2 ) + TS,(-£) n+1 (H (3n+1) / n -Y (3n+1) / n )/H n 

2 df m U ’ 31) 

5 m I- * “ = T h - \ + r l S 2 ( - f > n+1/Il2n ( 4 ’ 3 2) 


where 

J ,(x) -iru(-f) n /H 2n+1 (4.33) 

— r 

The final equations (4*22.) # ( 4 . 30 ) # (4*26) and (4*3l) are 
coupled through in and have to be solved simultaneously under 
the following boundary conditions I 


p, = 0# T, = T # T = T at X = -00 
•^1 1 o o 

P 2 = 0 at X = 0 


(4.34) 
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Ths Reynolds equations (4*22/4*30) and the mean 
temperature equations (4* 28/4*32) contain an unknown 
which is determined using the above conditions and the 
process mentioned earlier* Subsequently pressure and 
chc mean temperature can be obtained by solving equations 
(4.22-4.28) and (4*30-4.32) simultaneously- 

It may be noted at this stage that \** 0 corresponds 
to the case when inertia and convection terms have been 
ignored. For this case Reynolds and energy equations for 
the two specified regions reduce to 


dpq 

dX~ 

= J f (X) 




(4*35) 

dp £ 
dX~ 

= J_ f (x) 


Ti = 

t k + rS 1 f n+1 (H (3n+1, / n -y (3n+1)/n )/H ( 
n 1 

2n+l) (n+l) ,/n) 
(4*36) 

T = 
1 2 

T„+rm, ( - f ) n ' ! - 1 (H ( 3n+1 J - /n -y ( 3 n+1 ) /n, ^ 

n 2 

( 2 n+l) ( n-:-l ) /n) 


— — r Hrl /TT 2n 


T 

m l 

= T h + r l m l f / H 

(4*37) 

T 

m 2 

- , r-\ n+1 ~,2n 
= T h + r-^C-f) ,/H 



Equations similar to the last three sets of equations (4*35)/ 
( 4 * 36 ) / (4*37) were also obtained in chapter III. These can 
be solved by the process outlined earlier giving pressure 
and temperature distributions for the case X- 0* 
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koad and Traction * The non-dimensional load Wy and the 
non-dimensional traction/drag T p are obtained from the 
following expressions • 


W, 


-X - ,.-n 
1 m- Xr 
r 1 


H 


2n+l 


dX 


o 

/ 



m„x(-f ) n 


(4-38) 


T 


P 


' X 1 m, f 

• X 


: n 


H 


2 n 


ax 



H 2n 


dX 


(4*39) 


4*3 RESULTS AND DISCUSSION T 

The numerical values of the dimensionless pressure p 
and the mean film temperature T m are obtained by solving 
equations (4* 22) / (4* 28) and ( 4 . 30) / ( 4 * 32) simultaneously 
for various values of the flow behaviour index n(*4^ n £ 1*15) 
and squeeze velocity q(~0.09£ q£ +0.09)* The other 
dimensionless parameters which are important in the present 
context are A( = PU 2 a)/ p r (= Zk0a^ 2R~^ an d E^C =™^) . 

The values of these parameters X/ P r and are chosen to 
be 26xlO~ 5 / 5 2 50 6x10 _1 and 202xlo“ 6 respectively* The other 
representative numerical values used in this analysis are the 
following * 

U = 450 cm sec - ^/ h Q = 7x10*"^ cm/ R = 3 cm/ a = l«6xlO~^dyne~^an 2 * 
All the parameters/ especially n/q and X have profound effects 
on various bearing characteristics* This has been elaborated 
through various graphs and tables. It may be noted at this 
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stage that for the case a •* 0 / it was not possible to plot 
curves for all n and q, because the differences were not 
discernible in the graphs* So araphs have been made for 
X 4 0 mainly and curves for only a few values of n and q 
in case of X -* 0 have been drawn (denoted by 

The present case may be considered to be a generalization 
of the previous chapter as well as the work done -by Pascovici 
[156]* Pascovici / s case may be obtained directly from this 
problem by letting X ■* 0/ n = 1 in Reynolds and energy equations/ 
and ci = 0 in the empirical relation for m given in ( 4 * 4 )* 

Pressure Distribution '• The lubricant pressure p which is a 
function of n and q has been computed numerically by solving 
Reynolds equations (4*22) and ( 4 . 30)* The results have been 
presented in fig* (4* 2)* For both the cases of X / 0 and X •* 0 
(as mentioned earlier)/ the qualitative behaviour of p for 
the different values of q (for fixed n) is identical* This 
is in conformity with the results of Dowson et al [l2l] 
and those obtained in the first and the third chapters* A 
similar trend for both the cases is displayed by the pressure 
curves when n is varied and q is held fixed* This kind of 
behaviour was observed by Safar and Shawki [ 63 J for a thrust 
bearing and by Bukholz [73] for a journal bearing* The point 
of maximum pressure for the case X ^ 0 is of special Interest 
here* Mathematically speaking/ the point of maximum pressure 
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for the present case is the point of intersection of the 

curves I n (x) and J-(X;. This has been shown in fig* (4*1)? 

Physically it is the point at which inertia and viscous 

forces balance each other* The point of maximum pressure 

for the case A £ 0 is different from the point X = -X^ 

(at which — - is zero) and lies before X = -X, » However# 
o y r 

for the case X-* 0 / these two coincide i*e* the pressure 
gradient and the velocity gradient vanish simultaneously. 

A comparison between the two cases can also be made 
from the same figure* It is seen from fig* (4*2) that for 
each fixed n and q # the pressure curves for A 0 lie 
above the corresponding curves for X ** 0* One may therefore 
infer that the inclusion of inertia term leads to an 
increased pressure* It may be further noted that the inertia 
affect is less on the Newtonian fluid (n=l) in comparison 
to than for the dilatant fluid (n > l) • Increased pressure 
due to the presence of inertia has also been reportea by 
Safer [166] and others {,116}* 

Temperature Distribution i The mean film temperature T m is 
computed from differential equations (4*28) and (4*32)* 

The results for T m are elaborated through figs* (4*3/ and 
( 4 . 4 ) for various values of n and q (for both the cases 
X 4 0 and X -* 0) * It is clear from fig* (4*3) that for a 
fixed n and q(>0)/ T m increases as the lubricant flows 
towards the pressure peak region. Subsequently# it decreases 
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till the pressure is zero (ambient). For a fixed n and 
q^O / the behaviour of is similar to that mentioned 
above/ except near the centre line of contact (X=o) where 
it increases though only marginally. For a fixed n/ $ 
increases as q decreases whereas for a fixed q/ 5 increases 
as n increases* A similar trend for temperature with respect 
to n and q has also been observed in the chapter I. 

A comparison between the mean temperature for both the 
cases has also been made in fig. ( 4 . 4 )* It is observed from 
the figure that the mean temperature near the inlet for the 
case A 4 0/ is less than that for the case A-* 0. However/ 
the trend is reversed in the pressure peak region. One may 
interpret these results from a physical view point- For the 
case A ^ 0/ the convection dominant region extends farther 
than that for the case A-»0. Thus the heat lost by convection 
for A - * 0 is smaller as compared to that for A / 0. This may 
cause the temperature rise there. However/ as the lubricant 
flows towards the pressure peak region/ the heat transported 
by conduction becomes more significant for the case A**0- 
Consequently/ T m decreases continuously right up to the 
pressure peak region* Thus for the case A ^ 0/ T is higher 
than that for A**0 near the pressure peak region. 

Crmgistgncy Variation * The dimensionless consistency 5/ given 
by the empirical relation viz* 

p-tL+ta 

m s= m Q e^ m o 
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is depicted in fig. (4*5; for various values of n and q. 

It is evident that the consistency variation depends on 
the combined effects of pressure p and the mean film 
temperature T m « Since m for the pseudoplastic fluids (n< l) 
does not vary significantly with respect to q. Therefore/ 
for this case/ in is plotted for only one value of q(=-0.09). 
Moreover/ the consistency variation for all the values of 
n and q have similar trends i»e* in decreases near the inlet 
marginally because increase in T m is marginally higher than 
the corresponding increase in p* Thereafter/ in follows the 
pressure trend throughout the region i*e» in increases and 
then decreases near X = 0. For each q/ m increases with n 
throughout the region* However for each n/fri increases with 
q in the inlet region whereas the trend is reversed in whole 
of the pressure paah region. It is seen from the fig- (4*5) 
that variation in in for n ^ 1 (particularly n <l) is very 
less. It may therefore be noted that in may be treated as 
constant for pseudoplas oic fluids (n< l). 

The consistency variation for X -* 0 has also been 
presented along with that for X 0. The qualitative behaviou 
of in versus X (for all n and q) in this case is similar to 
that for X ^ 0. However/ quantitatively / m for X-* O is 
less than that for X /- 0 except in the pressure peah region 
where the trend is reversed. This is because of abnormal 
change in and p in that region which has already been 


mentioned 
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L oa d and Tractio n ' The surface force Wy and the traction 
force Tj, are calculated from the expressions (4*38) and 
(4*39) for various values of n and q and for both the 
cases X ^ O and \-*0* These results are presented in 
table (4*2) ♦ It is seen from the table that for a fixed q/ 

VI. and T„ both increase as n increases whereas for a fixed 
n/ both of them increase as q decreases* This is in accordance 
with the view point that whan the surfaces are approaching 
towards each other (i*e* q decreases) higher pressure is 
developed (see fig* (4*2))* Consequently/ the load and the 
traction both will be more* 

In order to study the effect of inertia on load and 
traction/ results for K~* 0 have also been presented in 
table (4*2)» It is seen chat the inclusion of the inertia 
force leads to increase in the load and the traction for 
all n and q. However/ the increase in results is not very- 
significant for the pseudoplastic fluids (n< l)* 

Sepa rate d R egions ) The Reynolds and the energy equations 

(4*22) / (4*28) / (4*30> and (4*32J contain (X a -X^ is the 

point at which vanishes)* This point X = -X^ separates 

3V 

the region £X £0 into two subregions -"iX £ ->g_ 

and 0* The values of for different n and q/ 

and for both the cases have been presented in table (4*1 )• It 
can be seen from the table that for a fixed n/ X^ increases 
as q increases whereas for a fixed q / X^ decreases as n 
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Increases* A similar trend was also obtained by Sinha 
et al [66] for X •* 0 and m constant- 

4-4 CONCLUSION * A theoretical analysis of cylindrical 
rollers moving with equal velocity and lubricated with a 
power law fluid has been presented- Thermal and inertia 
effects have also been taken into account together with the 
squeezing motion. A simple model for m given in equation (4-4) 
has baen considered. Average/mean values of the fluid inertia 
and the energy convection considered in the mometum and 
energy equations result in an anlytical solution for temperature 
and semi -analytical solutions for pressure and the mean 
temperature- 

Fluid inertia causes a significant increase in pressure 
and thus the load and the traction. This effect has even more 
significant in case of squeezing motion- This effect however/ 
does not have significant effect on the mean temperature. 

The fluid inertia tends to stretch the fluid film which 
results in a lubricatincr flow field larger than that of the 
inertialess case (X~*0) i-e- the point at infinity is shifted 
further away from the centre line of contact O- This also 
tends to shift the position "of pressure peak away from the 


origin 0 
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Table 4 

•1 • Values 

of x i 




q 


“X 

n / m o 

-0-09 

o.o 

*t 0 • 09 

A 

1.15/0.56 

0.507454 

0.532137 

0.564312 


1.00/0.75 

0.542656 

0# 564640 

0.593375 


0.545/86 

0.615625 

0.631562 

0.652031 

2xlO“ 5 

0.40/128 

0.656250 

0 .663531 

0,680562 


1.15/0.56 

0.528906 

0.553437 

0.585844 


1-00/0-75 

0.548830 

0.570569 

0.598645 

2 6x10 " 5 

0,545/86 

0.613262 

0. 634101 

0.654375 

0.40/123 

0,656731 

0 • 664062 

0. 681250 
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Table 4*2 . Normal force component and traction 


n / m o 

Wy 

T 

F 

q 



0.840294 

1.939247 

1 

O 

* 

O 

i 


1.15/0.56 

0.75196 7 

1.796557 

0.00 



0.661339 

1.648009 

4“0 * 0 9 



0.156338 

0.398510 

-0.09 


1.00/0.75 

0.142047 

0.372889 

0.00 

26xlO -5 


0.127634 

0.347025 

4-0 * 09 



0.053363 

0.176291 

-0.09 


0.545/86 

0-050922 

0.170916 

0.00 



0.048347 

0.165351 

4*0 * 0 9 



0.012929 

0.047981 

-0-09 


0.40/128 

0.012440 

0.046617 

0.00 



0.011919 

0.045219 

+0.09 



0.700509 

1.573574 

-0.09 


1.15/0-56 

0.632719 

1.481634 

0.00 



0.563006 

1.383317 

+0.09 


> 

0.148747 

0.374285 

-0.09 


1.00/0.75 

0.135689 

0.352230 

0.00 

2x10" 5 


0.122407 

0.329620 

+0.09 



0.051965 

0.169671 

-0.09 


0.545/86 

0.049722 

0.165095 

0.00 



0.047331 

0.160291 

+0.09 



0.012412 

0.044664 

-0.09 


0.40/128 

0.012027 

0.043809 

0.00 


0.011607 

0.042888 

+0.09 






Fig. 4.3 The mean temperature distribution T m against X. 



Fig. 4.4 Temperature distribution T m against X. 



CHAPTER V 


THERMAL EFPECTS IN LUBRICATION OP 
ROLLING/SLIDING CONTACTS 

5 . 1 INTRODUCTION *. 

Most of the roller bearing problems (pure rolling 

motion of cylinders) considering a power law lubricant 

have been studied using symmetric conditions i*e* both the 

rollers have the same dimension and velocity [ 30 / 33 / 57 / 66 / J* 

In all the preceding chapters also / the symmetric conditi 01 ^ 

for pure rolling were made* Particularly/ the lubricant fl° w 

about x-axis was made symmetric (see chapter l). How fat arG 

v 1 * o ^ ^ O 3^ 

these idealized symmetric conditions true under actual opa* c 
is in fact a matter of investigation. An obvious generali za ^° 
of the work presented earlier in this thesis would be the 
relaxation of the symmetric conditions. i.e* both the roll ers 
may be assumed to have different dimensions and velociti eS 
( anti- symmetric case)* Thus rolling/sliding operations come 
into picture* 

Several investigations on rolling/sliding concentrate^ 
contacts have been reported. Some of the early important vj°~^ s 
in this direction hava already been mentioned in the pre va -° us 
chapters [62/86/89/91/92 j * More recently Sadeghi and Dow 
presented a two dimensional solution to the problem of thermal 
EHL of rolling/sliding contacts* The Reynolds equation/ 
die two dimensional energy equation and the elastic equntr° n 
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have been solved simultaneously considering a pressure 
and temperature dependent viscosity for a synthetic 
hydrocarbon lubricant (XKM-109F) . Sadeghi et al [175] 
subsequently studied the same problem using an approximate 
method for the prediction of mid-film temperature and sliding 
traction* Rashid and Seireg [176] elaborated heat pa’tit .'o' 
and transient temperature distribution i n dry and lubr-i Cated 
rolling/sliding contacts- They Particularly commuted on 
the existence of the non-Newtonian behaviour of +he lub • -t 
in tho concentrated contacts- They fait that tvn „ v, 

tnis happens 

because the lubricant is subjected to extremal v -h-ir-rV, 

i nign pressures 

and shear-stresses which act for a very short time/ thus th- 
Newtonian property of the lubricant may not be valid. Rashid 
and Seireg [177] further mentioned that the severe thermal 
environment occuring in this type of contacts due to th^ 
combination of high pressure and sliding speeds is one of tho 
main factors in the malfunctioning of the machine dements 
such as gears/ bearings/ cams and drives- Bush et al [l78] 
analysed surface roughness effects in a rolli ng / slidi 1± 
contact using non-Newtonian lubricants such 

as a power law/ a 

non-linear sinh-law and a linear visco-dastic Maxwdl model. 
Expressions for bearing traction along with numerical resul ^ s 

have been presented- Apart from these/ there are s^vor 1 
other publications on lubricated rolU ng/sliding Uno contacts 

where Newtonian and different non-Newtonian models hav^ b^n 

considered in the analysis- Pressure and temner^rnr-. * 

^emperarure depondenc 
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on non-Newtonian and Newtonian viscosities of the lubricants 
have also been taken into account* 

In most of the references cited above/ anti-symmetric 
conditions have been used* It is also clear that the lubricant 
properties change with pressure and temperature during rolling/ 
sliding concentrated contacts* Particularly/ the lubricant- 
viscosity changes significantly. The non-Newtonian character- 
istic of the lubricant also seems to be important. Since 
the lubricant pressure and temperature in case of rolling/ 
sliding are expected to be higher than those for the pure 
rolling case (as considered in the previous chapters)/ it 
becomes far important to. incorporate the non-Newtonian behaviour 
of tho lubricant together with variable viscosity/consistency. 

In the present analysis/ a power law fluid film 
lubrication of cylindrical rollers having different dimensions 
and velocities/ Is studied under cavitation boundary conditions* 
Since the rolling bearing geometry is equivalent to a cylinder- 
plane geometry/ the latter is prefered for mathematical 
calculation. Lubricant consistency is allowed to vary in a 
manner similar to that' mentioned in chapter III* The lubricant 

I 

inertia which is generally small has been ignored / as this 
is likely to make the solutions extremely complicated* A very 
simple approach is adopted in the present analysis in order 
to separate the zones for positive and negative velocity 
gradients* 
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5.2 MATHEMATICAL ANALYSIS * 

The governing equations with the usual assumptions of 
hydrodynamic lubrication for the problem under consideration/ 
in which the two rollers are assumed to have different 
dimensions and velocities/ are 

an) (5.1) 

dx 3 y 3 y 3 y 



(5.2) 

(5.3) 

(5.4) 


From the point of view of mathematical analysis/ the contact 

between two cvlinders of radii R and R. / as shown in fig- 

a .0 

(5 -1(a)) can be described by an equivalent cylinder near 
a plane/ shown in fig. (5.1(b)) [126]. 


Film thickness h for the geometry given in fig. (5.1(b)) 


is expressed as 

h =h o 

where 



(5-5) 


1 11 

R ~ R + R, 
3. O 


(5-6) 
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The boundary conditions for the governing equations (5.1) 
and (5.2) are as follows * 


u 

to -1 

11 

9 V t= 

0 

at y 

= 0 

(5. 

.7) 

u 

= U 2 

/ V = 

U 2 

dh 

dx 

at y = h 

( 5 ‘ 

» 8) 

p 

= 0 

at x = 

= — 

OO 










(5. 

9) 

p 

= 0 / 

aS = 

0 

PJ 

C+ 

X 

= x 2 




In order to solve the' system of equations (5.1) and (5.2) 
for various bearing characterics/ one has to first of all 
determine the sign of the velocitv gradient i.e. to 

ay 

eastablish the regions of positive or negative velocity 
gradients. Approximately/ these regions can be obtained 
easily for the Newtonian fluid (n=l) and it is expected that 
there will not be a significant deviation for non-Newtonian 
fluids* 

Equations (5.1) and (5. 2) can be solved for n = 1 and m to be 
a constant or a function of x-alone / one can obtain [89] 


u - A + (U 2-V £ - (h -y> 

= 6(u 1 +u 2 )m(h-h L ) /h 3 


(5.10) 

(5.11) 


It can be seen that u is a linear function of y at points 
where = 0 . At these points/ the velocity gradient 


a* (u 2" u i } 

ay - K 

It can also be seen from fig. (5.2) that for each x/ can 

BY 


/ which can never be equal to zero/ since U 2 ■£ 
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vanish at only one point of y(=6) in the regions -«> ^X< -X^ 

and -*x^ < x< x^* Thus these regions can be divided into 
four sub-regions separated by 5-profile/ as shown in fig*(5»2)» 
Let ^ /U 2 an< ^ U 3 /U 4 " i ' ne va l oc i“2- es the respective 
regions (see fig. (5. 2))/ then according to the velocity 
profile/ shown in the figure / one may Observe the following ; 

/-» sE X < (5.12) 


3Y 

3u, 

3Y 


a u. 


V 


S o / 5 ^ y ^ h 


^ 0/ 0 ^ y g 5 


3 u, 

d~Y~ 


0 at v = 6 


3 u 3 

aT 

9_^4 

ay 


a u 


Z 0 / 0 £ y£ 6 


gO, «S yjh 


_4 _ 3 5s 

ay ay 


= 0 at y = 5 


i-xx< x < x 2 


(5.1?) 


It is expected that the nature of the velocity profiles 
for different n as well as variable m may remain unchanged 
but 6 will certainly differ in each case- 


Using the sign for the velocity gradients given in (5.12) 
and integrating equation (5.1) for the region 5 g y g h/ one 
can obtain 


3 u^ 
dY~ 


= ( 


^,Vn (y ^,Vn 
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Further integration of this, using condition ( 5 » 8 )/ gives 

dp 


{(h-6) (n+l) ' /n ~(y-5) (n - f - 1)/n J 

1 


U-, = U„ - TrUr“ srr-) ' Hh-d ; ' ~-Ky- 

(5.14) 

Similarly for the other regions/ the following expressions 
can be obtained i 

l l 

dp 

‘2 

dp 

l 2 


u 2 = U]L _(J2j->(±- a n.) 1 /n f6 (n + l)/^ u . y) (n+l)/ n , (5 . 15 

h-C^H- h S i ) 1/n !S <n+1,/n -<5-y) (n+1)/n ! Cs.is) 


u 3 ■ U l" rv n+I 




(5.17) 


Now the volume flux Q for the region (-« g x < -x^) is 
obtained as 


Q = / u dy 
o 


or 


6 h 

Q = / Uj dy + / Ui dy 

° " 5 


>r Q = u 1 6 +U 2 (h- 6 )-( 2 |^) (~ { 6 (2n+l) / n + (h- 6 ) 


( 2 n+l)/n 
(5.18) 


Since the flux Q is constant through the every cross section/ 
thus the flux Q may be equated to the flux through the point 
x = ~x^ [179]. This point is selected in order to enforce 
the continuity of flow. 
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The flux Q through the point x = ~x^( where ^ = 0) can 
be evaluated directly from the equation (5.1)/ yielding 

(U, +U_) 

Q ( -x^ ) = h x (5-19) 

Finally, equating the fluxes given by (5.18) and (5-19) and 
simplifying/ one can obtain the following Reynolds equation 


dpj, 

asr 


= m. 


U +U 

2n+l> n r (U 1 -U 2 )5+ u 2 h - <-4 r J,h 


/ \ n r 

. n ; 


^ ( n+1 ) /n ^ h—6 ) ^ ^ ^ 


- — in - - ~ £ x ^ - 


*1 


(5.20) 

Similarly for the other region/ the Reynolds equation may 
be written as 


U, + U 0 

dp, 2nt . n ( “2-- ,h l- U 2 h - (U l' U 2 )6 „ 

cBT = " m 2 ( “n : “) t rSntTTTnTT^ TTTIn+TTTn^ 7 " X 1 - xS x 2 


+ ( h-6 ) 


Use of the matching conditions: 


U 1 = U 2 


at y 


(5.21) 


(5. 22) 


u. 


u. 


into (5.14-5.17) yields 

dp 


U 2"’ U 1 " ( n +r" )( rr (h-6) (n+l) ' /n -6 (n+l)/n i= 0 / 


"3x 


•co ^ X < -Xj_ (5.23) 

U 2“ U 1 +( H&:" ) (“ MJ d^" )1/n ^ n+1 Vn_ 6 (iM- 1 )/nj s _ ° # 

-X, < X < X, 


(5.24) 
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dp l dp 2 

Eliminating and from equations (5*23) and (5-24)/ 

using Reynolds equations (5-20) and (5. 21)/ one can obtain 
a single relation as follows • 


U, +U, 




{ ( h-S ) ( n+1 } ' /n -6 ( n+1 ^ ' /n l { (Uj -U 2 ) 6 +U 2 h- (“ 2 “) hj} 


[5 (2n+l)/n + (h ^ 6) (2n+l)/nj 


- C 


(5. 25) 


It must be noted here that the equation (5.25) Cannot be 
used to evaluate $ at x = -x^ and x = x 2 since £ 0 


there* 


Energy equation * The energy equation for this problem is 
considered to be the same as that considered in chapter III i»e« 


3 2 T _ f rru 

“T 


t^H - 1 

3 y 


2 


ay~ " 3Y 

where the heat of convection has been neglected* 


(5-26) 


The boundary conditions for this equation are 
T = T q1 at y = 0 
T = T h at y = h 


(5,27) 


Integration of equation (5.26) twice for the region / 
- ^ x ^ -x-l / gives 


T 11 = 


n 


(j “) ("“ gj~) ( n+1 } /n (y-6 ) ( 3n+1 ) /n +c 1 y+d 1 

(5.28) 


T2n+lJT3n+l) 


n 


■12 ~ “ TTn+lTT 3 n+lT v h 


ra l^ .I dp ^ (n+1) /n/ e \ (3n+l)/n 


(^) ( 


clx 


(5-y) 


+ c 2 y+d 2 
(5. 29) 
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where T 11 and are the lubricant temperature in the 
regions/ 6 £ y ^ h and y £ 6 respectively/ c-^/d^ and 
c 2 ,< ^2 are the intsura'kion constants* 

Use of the temperature matching condition 

T n = t 12 at y = 6 

and the matching heat flux 


3T 

3Y 


11 


•- = k 


3T 

-g-y 


12 


y = 5 


implies 


c^ = c^ = c(say) 

= d 2 = d( say) • 

The constants c and d are evaluated using the boundary 
conditions (5*27) and are given below • 


T, -T 


= -TT— + ‘Tar’ <FT tn+1)/n [(h-6) <3n+1) /n. 6 <3n + l) / 


dp. 


m, a_ 


d = T 


“1 

dp 


(5.?0) 


ol 


+ (JJUS) ( L_ ^Li) ( n+ l ) /n 5 ( 3n+l ) /n 


dx 1 


where 


A n = 


n* 


1"2n+l) 1 3 n+1 ) 


Thus T^^ and T^ 2 are explicitelv known functions of x and 

Y* giving temperature distribution in the region under 
consideration. 
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Finally the mean temperature T m t as defined in equation (5*4)/ 


may be written as 
1 h 

T = r / T dv 
m l h o 


T ni = h c i T i2 d y - / T n J 


T,+T , S, (3n+l)/n (3n+l)/n 

orT mi = Jl^.21 + [(h-6) 


+ 6 


] 


nS. 


( 4 n+1 ) /n (4 n+1 ) /n 


- r4^+iTh { (h " 6) +5 } 


(5-31/ 


where 


, m l A n^ ,1_ dp J^ ( n+1) /n 
S 1 " ( k } ( -m 1 dx ' 


(5-32) 


Proceeding in the same way / the mean temperature for the 
region -x^ ^ x £ x ^ t can also be written as 


T„ = ( 


m 


T, +T . S 0 (3n+l)/n 

-S °i) + -£[(1,-6) +5 


(3 n+1) /n 

' ] 


nS, 


( 4 n+1 ) /n ( 4n-i-l ) /n_ 


- T4^?lTh^ (h " 6) J 


(5.33) 


wnere 


m A . dp (n+l)/n 

i 2 rx\ t l 2-, 

- C- x ~)(- — ^) 


(5.34) 


Using the non-dimensional scheme given in the previous 
chapter together with the following • 

T m = f!T m - 6 = 6/h c , 0 = U x /U 2 (U x > U 2 ) (5.35) 

equations (5*20)# (5«2l)/ (5*25)/ (5«3l)/ (5.33) may 
be rewritten as 


dp. 


■1 
dX" 


n 


(5*36) 
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d Po _ . - , 

wr = ~ m 2 -f x 


■J. n n 




(5.37) 

(n+l)/n (n+l)/n ,l+tJ\ i 

(h-J) -6 HW-DS-m-firH! 


= 0 


|- 6 (2n+l)/n + (H ^) <2n+l)/nj 


T +T , 
n ol 


2 


+ S 1 d n 7 s ( V 


n+1 - 




1 ? *1 714-1 —* 
5 m 2 * - S 2” + 5 2 d n V ' s ( "V a ' 


(5« 38) 
(5.39) 

(5*40) 


where 


Po~T m _+T. 


m l = m Q e 


Pl” T m 1 +T o - - p 2 m 

s 1 / m 2 = m Q e 


, U 2 h o p w h oa/2 

r = 


d n = n/((3n+l)2 n+1 ) 


5 = r 

X L 


(U-l)5 


-] 


(5-41) 

(5*42) 

(5*43) 

(5*44) 


*X 


£-(2n+l)/n + ( H _£-) (2n+l)/n 

(3n+l)/n (3n+l)/n 

= |{(H^6) +6 1 

{(H-ff) (4n+l) /n + £ (4n+l) /n ] 


n 


(4n+llH 


Equations (5.36), (5.39) and (5*3?) , (5-40) have to be 
simultaneously using the following conditions * 


p l 


= 0 at X = -« 


(5-45) 

solved 

(5.46) 
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p 2 = 0 

dp, 

aor = 0 


at X = X 2 


(5.47) 


From equation (5.37) * it may be observed that the 
condition 
dp, 

wr =■ 0 at x = 

implies 


(U-l) 6 2 + h 2 


(u+i) tt 
“ H 1 


= O 


(5*48) 


Since 6*2 comes out to be H 2 /2 (for detail* see numerical 
calculation* section (5«3)) * it follows that 

X 2 = X x (5.49) 

Thus X = X-^ gives the point of cavitation. 

5.3 ‘ RESULTS AND DISCUSSION I 

Numerical solutions of the Reynolds equations (5.36*5.37) 
and the energy equations (5.39*5.40) have been obtained 
for anti -symmetric * viscous incompressible flow of a power 
law fluid through the gap between a cylinder and a plane 
(see fig* (5.1)). The results of this investigation are 
assessed in terms of parameters n and U(= U^/U 2 * U^> U 2 ) * 

The values of n have already been mentioned in earlier chapter 
The sliding parameter U is chosen to lie between 1 and 
1.2 (i« 3 * max of 20% slip)* The parameter U arises due to 



12 ? 


the consideration of anti-symmetric conditions and is 
important because the presence of sliding (U> l) is likely 
to yield more pressure and temperature as compared to that 
for pure rolling (U=l). The significance of U along with 
n has been demonstrated, through table and graphs* For the 
numerical calculation/ the following representative values 
have been used 

U 2 - 400 cm sec"' 1 '/ h = 5x10*"^ cm/ oc - 1.6xl0~ 9 dyne**' 1 ' cm' 1 '/ 

R = 3 cm/ T^ = 1*5/ T q ^ =1* The thermal factor 
f s = (t^h^p ^fc 0 /2R/ka) is chosen to be 5* 

An important feature of this problem is that this reduces 
to (a) the exact problem done in [ 33 ] for and 

m = constant / (b) a similar problem done in [l5o] for 

n = 1 and a - 0 and (c) the exact problem done in chapter III 
for Uj = U 9 . 

Nume r ical Solutions ♦ The Reynolds and energy equations 

are coupled through m and contain two unknowns 5 (the locus 

of points at which = 0) and X, (X=-X. is the point of 

0 2 «*■ 

maximum pressure) * These unknowns are also present in 
equation (5*38)* As there is no symmetry (U ^ l) 
for 6 -distribution / it is necessary to solve equation (5*38) 
for § over the complete region under consideration# along 
with Reynolds and energy equations* The actual process 
followed for the numerical calculations is given below* 



The algebraic equation (5.33) in <5 contains X and 
explicitly- First of all an initial value of X is 
assigned i-e. the point a~ minus infinity is replaced 
by a large but a finite negative value- An arbitrary value 
of X^ is chosen- The value of 5 at the inlet is obtained 
by solving the algebraic^ equation (5-38) for 6 using the 
bisection method with a reasonable tolerance (say 10 ) - 

These values of X and 5 are substituted in the energy equation 
(5-39). T m at the inle-c is then obtained from solving the 
algebraic energy equation by prescribing m 0 and p^ (=0 at 
X =■: ) and using the earlier numerical technique- The 

sa^s X/ m Q and the computed 5 and T^ are used in the 
differential equation (5-36)* Fourth order Runga-Kutta 
method is used to evaluate p^ at X = X+AX - For this value 
of X (=X+AX) / equation (5.38) is solved for 6- These new X / 

6 and p-^ are substituted in the energy equation which again 
yields T m , X and § are used to calculate p^ as a solution 
of equation (5-36) at another point X = X+2AX. This process 
is repeated for further values of X till such time that 

l 

'equation (5-38) yields values of 6 satisfying 0 E ^ H 
(because 5 cannot exceed H) - It may be noted that at X = -X^ / 
<5 does not exist (see previous section). Further in the 
neighbour-hood (n-hood). of X = -X^ / the determined values 
of 5 do not lie in the interval 0 £ 5 £ H- Hence in the 
n hood of X = — X^ / 6 (=say 6 ) has to be determined solely 
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on the basis of physical consideration* It should be 

_ * 

emphasized here that 5 does not refer to the locus of 
points at ■which the velocity gradient vanishes* Let one 
consider an e -n # hood (-X^-e ^ X £ -X^+e/ e>0) of the 
point X = -X x where there exists no <5 lying in the interval 
0<j 6 ^ H and satisfying -equation (5*38). To ease the 

mathematical complexity/ a .linear profile for E given 
below 

6 *(x) = - ~ (X + X ± -e) (5*50) 

is assumed • This profile was chosen in such a manner that 

it not only satisfies the Reynolds equations at X = -X^ / 

but also makes the pressure curve continuous* Having 

determined <5 = E* using equation (5*50) in the n'hood 

of the point X = -X^ / the procedure outlined earlier is 

adopted to evaluate p^ and 

and p„ / 5 for the reaion -X, ^ X £ -X, +e. Since in the 
2 m 2 "1 1 

interval -X^+e ^ X ^ X 2” e * $ (roots of equation (5*38)) 

satisfies the inequality OS 5 ^ H/ so E in the region 

-X 1 +e ^ X^ X ? -e is determined (by the same procedure followed 

near the inlet) along with p 2 and T m * Subsequently in the 

interval X 2 -s ^ X ^ X 2 (since equation (5*38) is not valid 

at X = xJ / the same type of linear profile given below 
* H, 

5 (x) = ^|(x-x 2 +e) / x 2 -e£x£x 2 

is assumed (The complete 5 -profile has been shown in 
fig. (5-2))* 


T 


m. 


for the reaion 


-X, X£ -X, 


-L 


(5.51) 
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This will satisfy the Reynolds cavitation boundary 

dp 

condition (5*47) i.e* = 0 at X = (5.52) 

iff X 2 = X ± 

5 in the region X^-s ^ is calculated using 

relation (5*48) together with. p 2 and T m * If the computed 
value of p 2 at X = X^ satisfies the condition (5.47) 

(i*e. p 2 =0 at X = X^) / the assumed arbitrary value of 

X x was correct* Otherwise another value of X^ is assigned 
and the whole process is repeated so long as p 2 vanishes 
at X = X^. Thus X^ and 6 (as a function of X) are computed 
along with pressures Pj'^2 sn ^ ’ fche mean temperatures T /T^ 
The complete pressure and the mean temperature profiles 
as functions of n and U have been presented in the following 
sections ♦ 

Pressure Profile * The pressure profile p has been presented 
for various values of n and U in fig. (5.3). It can be seen 
that for each fixed n and U/ the pattern of the pressure 
profile with respect to X is exactly similar to that shown 
in previous chapters* It can also be noted that for a fixed 
U# the qualitative behaviour of p with respect to n remains 
unaltered (see earlier chapters) i»e* p increases with n. For 
a fixed n. p increases with U* This can be seen through the 
fig. (5.3) for values of n £ 1. p depicts the same trend 
for n< 1/ also/ however it has not been shown because of 
overlapping of curves* p increases with U which implies that 
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hydrodynamic pressure due to sliding is higher than that 
for the pure rolling* This is justified also because the 
rolling friction is always less than the sliding friction- 

T emperature Profile * The mean temperature distribution 
has been depicted in fig* (5-4) as a function of n and ti- 
lt is clear that the qualitative behaviour of the mean 
temperature T m versus X is quite similar to that obtained 
in chapter III* For a fixed U / T_ increases with increasing 
value of n* This pattern also is in accordance with the 
results mentioned in that chapter (for q = 0) * For a fixed 
ru T m increases with U (its behaviour is similar to that of 
the pressure profile! • Hence the same geometrical interpretat 
ion may follow here* It may be easily concluded that the 
sliding temperature is higher than that for the pure 
rolling- 

Con sis tency Profile • The characteristic behaviour 
of the consistency in i»e* over all change in in due to the 
combined effect of pressure p and the temperature T m (see 
relation (5*41)) has been demonstrated through the fig* (5*5) ■ 
It is seen that the consistency distribution in versus X is 
qualitatively almost similar to the pressure profile* The 
difference appears (in case of n > l) near the inlet as well 
as outlet* However / for n£ 1/ the difference appears only ; 
near the outlet* Near the inlet (for n > l)# in decreases 
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because the increase in p is unable to compensate for 

the temperature rise* After that it is almost flattened 

because of equal contribution of p and T m * Again at the 

outlet (for all n) / m increases (though very marginal) 

since the impact of the reduced temperature is higher than 

that of the reduced pressure? Further the consistency 

peak in the low pressure gradient region is quite sharp* 

This happens due to sharp increase in p and decrease in 

T m (when X -» X^) simultaneously* For a fixed U/ rn increases 

with n* Similarly for a fixed n(> l) / it increases with ti 

except near the inlet where the trend is reversed* For a 

fixed n(£ 1)#- a similar trend has been obtained/ however 

it has not bo shown due to overlapping of curves* 

♦ 

kpa d, a nd Traction • Load Wy and traction Tp have been computed 
for various 'values of n and U and have been presented in 
table (5*1). For a fixed value of U/ Wy as well as T p both 
have similar trend with respect to n as discussed in the 
previous chapters* For each fixed ru Wy and T p increase 
with U. Hence it follows that load £md traction increase 
due to sliding. 

^P- l . m Rupture - Film cavitation points have also been calculat 
and have been displayed in the same table (5*1)* It may be 
seen from the table that for a fixed value of U / the cavitati 
points follow the same trend as mentioned earlier (i*e* X^ 
increases as n decreases)* Similarly for each fixed n* X^ 
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increases as U decreases* It seems that cavitation during 
the sliding motion occurs earlier to that of the pure rolling 
case* 

5.4 CONCLUSION : 

With the usual assumptions of hydro dynamic lubrication/ 
an endeavour is made to study thermal lubrication of rolling/ 
sliding contacts by an incompressible power law fluid. The 
thermal Reynolds and energy equations which are functions 
of consistency (m) / sliding velocity (u) and the location 
of points at which the velocity gradients vanish ( 5 )/ are 
derived. Semi-analytical solutions for pressure p and the 
mean temperature T m are obtained for the isothermal boundaries* 
It may be concluded from the results obtained that 

(a) there is a significant increase in pressure (hence 
the load and traction) with flow index n for a fixed value of 
U. A similar' trend follows when U is varied (from 1 to 1.2) 
and n is h£»ld fixed. 

(b) there is also a significant change in the mean temperature 
with respect to n and U* 

(c) In addition / the consistency variation resulting from 
p and T m is also important- Particularly/ this variation 

is quite significant in the pressure peah region for Newtonian 
and dilatant fluids* 



Table- 5-1 * Bearing characteristics for various values of 
n and U 


n/rn 
' o 

U 

X 1 = X 2 

W 

w y 

T 

F 


1-20 

0.420500' 

1.802456 

1.160597 


1.15 

0-423172 

1.789505 

1.148427 

1 * 15/0*5 6 

1.10 

0.425062 

1.774255 

1.133527 


1.05 

0.4267 37 

1.756257 

1.115375 


1.01 

0.427515 

1.740753 

1.100124 


1. 20 

0.457812 

0.460492 

0. 253302 


1.15 

0.461625 

0.453943 

0.250487 

1.00/0.75 

1.10 

0.463750 

0,446638 

0, 246039 


1.05 

0.464625 

0.439343 

0. 241624 


1.01 

0.465000' 

0.433326 

0. 237 861 


1-20 

0.474590 

0.157524 

0.076426 


1.15 

0.481500 

0.156369 

0.07 6375 

0 • 545 /86 

1.10 

0.4857 25 

0 . 1^5055 

0.076145 


1.05 

0.489300 

0. 153657 

0.075812 


1.01 

0.495125 

0.152522 

0.075730 


1.20 

0.483750 

0.037 37 2 

0.0175410 


1.15 

0.485000 

0.037140 

0.0175131 

0.400/128 

✓ 

1.10 

0.495750 

0.036910 

0.0174810 


1.05 

0.511500 

0.03667 2 

0.0174570 


1.01 

0.514062 

0 .036458 

0.0174290 




Lubrication of two cylindrical Fig. 5.1(b) Lubrication of a cylinder 
rollers having different dimensions on a - plane 

and velocities. 
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Fig. 5.5 The consistency distribution m vs.X 
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